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• Abstract 

We investigate the basic features of the gluon density predicted by a renormaUsation 
. group improved small- a; equation which incorporates both the gluon splitting function at 

' leading coUinear level and the exact BFKL kernel at next-to-leading level. We provide 

resummed results for the Green's function and its hard Pomeron exponent a's(as), and 
for the splitting function and its critical exponent LodcXs)- We find that non- linear resum- 
mation effects considerably extend the validity of the hard Pomeron regime by decreasing 
diffusion corrections to the Green's function exponent and by slowing down the drift 
towards the non-perturbative Pomeron regime. As in previous analyses, the resummed 
exponents are reduced to phenomenologically interesting values. Furthermore, signifi- 
cant preasymptotic effects are observed. In particular, the resummed splitting function 
departs from the DGLAP result in the moderate small- cc region, showing a shallow dip 
followed by the expected power increase in the very small-x region. Finally, we outline 
the extension of the resummation procedure to include the photon impact factors. 



1 Introduction 



Progress in understanding small-x physics has been characterized by quite a number of steps: 
first the BFKL evolution equation and its early prediction of the small-x rise of hard cross 
sections, leading to the notion of hard Pomeron in perturbative QCD; then, the qualitative 
confirmation of such a rise at HERA showing however a somewhat milder effect and, 
at the same time, good agreement with DGLAP evolution [S] at two-loop level; then the 
parallel calculation of the next-to-leading (NL) BFKL kernel [^[3]) leading to a dramatic 
decrease of the effect and to possible instabilities p l| 7 1 |5j of the leading log s series; finally, the 
proposal of various resummation approaches [HlEl^DElEl^SEl and recipes to stabilize 
the series, in order to provide reliable predictions for processes with two hard scales and 
DIS-type processes. 

The resummation approach proposed by some of us [Hl llUlITT] and summarized in Sec. |21 
identifies a few physical QCD effects that lead to large corrections. Firstly, the cross sec- 
tion dependence on the ratio of the hard scales of the problem, which is constrained by the 
renormalisation group (RG) requirement of single-logarithmic scaling violations in the rele- 
vant Bjorken variables. Secondly, the occurrence, at NL level, of the non-singular part (in 
moment space) of the anomalous dimension, yielding a sizable negative contribution. Finally, 
the running coupling effects which modify and make ambiguous the very notion of a hard 
Pomeron. 

A key effect of the running coupling is that the BFKL evolution drifts towards smaller 
momentum scales, which are more strongly coupled, thus making non-perturbative physics 
more important at high energies. This means that the asymptotically leading Pomeron 
LOp is actually a non-perturbative strong coupling quantity This feature can be 

taken into account by the initial condition in the DGLAP evolution of structure functions, 
but may be a problem in the processes with two hard scales (like Mueller-Navelet jets |19| . 
7*7* scattering 20^ etc.) where the perturbative hard Pomeron behavior can be observed at 
intermediate energies only. 

Recently, it has been noticed that the transition to the Pomeron regime is driven, in some 
small-x models, by a sudden tunneling effect [22122] at moderate values of as{t)logl/x, so 
that the 6-expansion [^S] may be needed to suppress the Pomeron and to identify the hard 
Pomeron exponent ujs{t) and its diffusion corrections |24 | l25 | l8ll26j (here t = logfc^/A^, where 
k is the transverse momentum of the hard probe, and A = Aqcd)- Furthermore, the gluon 
splitting function is expected to be power behaved in the small-x region too, but with a differ- 
ent exponent uJc{t), due to running coupling effects. Therefore, in a resummed approach with 
running coupling one has to investigate various high-energy exponents: the hard Pomeron 
index ujs{t) just mentioned, the resummed anomalous dimension singularity uJc{t), which are 
generally different and perturbatively calculable, finally the asymptotic Pomeron ujp which is 
determined by the strong coupling behavior of the model. 

The calculation of cug and lOc was performed in the renormalisation group improved (RGI) 
approach of jJJ. The result was that ^^(t) carries important non- linear effects, leading to 
a stable and sizable decrease with respect to its LL BFKL value, and that ujc{t) is sizably 
smaller than ujs(t) also. However, the method of solution of the RGI equation used in 
was best suited for the homogeneous equation, rather than the Green's function (cfr. Sec. [21) • 
Therefore, no real estimate of hard small-x cross sections was really possible. 

The purpose of the present paper is to further investigate the RGI approach by pro- 
viding a numerical calculation in k and rapidity space of the Green's function and of the 
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corresponding splitting function. By then using fc-factorization and the corresponding 
impact factors 28,29,30,31], this sets the ground for a full cross section calculation. Here 
we also provide the high energy exponents and a semi-analytical treatment of the diffusion 
corrections. Part of the results of this paper have been summarized elsewhere j32j . 

In order to perform such an analysis, we introduce a resummation scheme slightly different 
from that proposed in which turns out to be more convenient for numerical implemen- 
tation, and belongs to a class of schemes that are identical modulo NNLj; (and NLO in Q^) 
ambiguities intrinsic in the resummation approach. Recall, that — as summarized in the 
introductory Sec. |2l — the RGI approach incorporates leading and next-to-leading kernel 
information exactly, with some extra w-dependence {to is the Mellin variable conjugated to 
Y ~ log 1/x) so as to implement the RG constraints and the resummation of leading log 
collinear singularities mentioned before. Such requirements fix the form of the w-dependence 
of the kernel, apart from NNL terms, which remain and allow some freedom in the choice of 
the resummation scheme. 

The exact definition of the kernel and of the resummation scheme is provided in Sec. Ol 
Stated in words, the main difference of the present formulation with respect to that of Ref. jllj 
is that the resummation of the collinear behavior quoted before is obtained here by the uj- 
dependence of the leading kernel, rather than by a string of subleading ones. This allows us 
to include the full w-dependence of the one-loop anomalous dimension in a more direct way, 
while of course, leading plus NL kernel information is correctly incorporated, as in all such 
schemes. 

The detailed investigation of the gluon Green's function with its hard Pomeron behavior 
and its diffusion corrections is performed in Sec.0J by analytical and numerical methods. The 
full numerical evaluation relies on the method introduced in Ref. |33i . Through the numerical 
study we are able to analyze the border between perturbative and non-perturbative Pomeron 
behavior, at realistic values of Y and a^, and to extract the leading terms (~ bY, ~ b^Y^) 
in the exponent of the perturbative part. Such terms can also be calculated analytically 
by the 6-expansion method fMi 23 . We are thus able to identify both the hard Pomeron 
exponent at order 0(b) and its diffusion corrections, and we notice sizable non-linear effects 
which stabilize the intercept, decrease the diffusion effects and slow down the drift towards 
the non-perturbative Pomeron regime. 

We also provide in Sec. [5] the resummed splitting function. At the analytical level, we 
notice that the w-expansion method |inilllj allows one to define a resummed characteristic 
function which, in the saddle-point approximation, can be related to the "duality" approach 
of Ref. depending on the choice of the intercept in the latter. Beyond the saddle-point 
estimate, the resummed splitting function is evaluated numerically by the method of Ref. j34j . 
and shows a power increase ~ x~''^''^"''^ in the very small-x region, together with a shallow 
dip (compared to the DGLAP result) at moderately small x values. 

A preliminary discussion of the off-shell photon impact factors is provided in Sec. IHl Here 
we show how the resummation scheme incorporating collinear leading logs can be extended 
to the impact factor, and how the latter can be extracted from the result obtained in the 
recent literature We finally summarize and discuss our results in Sec. [3 
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2 Renormalisation Group improved approach 



The size of subleading corrections j^l^ to the BFKL kernel K,{k,k') and the ensuing insta- 
bihties jni|7||S] make it mandatory to understand the physical origin of the large terms and 
possibly resum them. In a series of papers ^9, .10,llj (for a review see ISSJ) it was argued that 
most of the large corrections were due to collinear contributions, so as to achieve consistency 
of high-energy factorization j^Zj at subleading level with the renormalisation group. This 
requires resummation ^ of both the energy scale-dependent terms of the kernel j5_^ and of 
the leading- log collinear logarithms ^10^ for both Q ^ Qo and Q <C Qo) with Q, Qq being 
the hard scales of the process. In the following we summarize the approach of which 
incorporates both the renormalisation group requirements and the known exact forms of the 
leading ^ and next-to-leading [^lE] BFKL kernel. A resummation for anomalous dimen- 
sions within a single collinear regime Q ^ Qo has been proposed in ^2]) and alternative 
resummations in [T^ITllIT^ . 

2.1 fc- factorization and high-energy exponents 

We consider a general process of scattering of two hard probes A and B with scales Q and 
Qo at high center-of-mass energy y^. We assume that the cross section can be written in the 
following fc-factorized form J.7<: 

where and are dimensionless impact factors which characterize the probes and ensure 
that \k\ (|fco|) is of order Q (Qo), and the gluon Green's function is defined by 

gUk,ko) = {k\[uj-IC^]-^\ko) . (2) 

The function /C^; is the kernel of the small-x equation of the general form 

ujgUk,ko) = 6\k-ko) + J ^ICUk,k')gUk',ko) . (3) 

The factorization formula involving two-(Regge)gluon exchange, has been justified up to 
NL logs level in Refs. 28. for initial partons and in j29u30) for physical probes. At further 
subleading levels, many (Regge)gluon Green's functions contribute to the cross section as 
well, due to the s-channel iteration. However, our purpose here is to incorporate leading- 
twist collinear behavior, and at that level the two-gluon contribution is dominant, so that we 
shall consider only the contribution in the following. 

While fc-factorization is supposed to be valid for as < w ^ 1, we shall sometimes extra- 
polate Eq. (P) to sizable values of w = 0(1) and moderate values of s, encouraged by the 
stability of our resummation, and by the possibility of incorporating phase space thresholds 
in Eq. (pQ) (cfr. Sec. E}. It should be kept in mind that such a region lies outside the validity 
range of Eq. (P), so that the extrapolated Green's function loses — most probably — its 
original meaning as two-(Regge)gluon propagator. 

In writing Eq. we have performed the choice of energy scale sq = QQo, in terms of 
which the high energy kinematics shows a simpler phase space, as explained in more detail 
in Sec. Actually, for intermediate subenergies it is more convenient to introduce as energy 
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variables the scalar products of type v = 2k^kQ, which have |fc||fco| as threshold, so that 
|fc||fco|/i^ is a good Mellin variable. Correspondingly, the energy dependence of the Green's 
function and of the impact factors is defined by [k = \k\,kQ = |fcol) 

g{,^,k,ko) = [p: fyi^y a.(fc,fco) (4) 



and 



2-111 ykko 



In this paper, we are mostly interested in the properties of the two-scale Green's function 
and of its high-energy exponents. It was pointed out in that, in the improved approach 
with running coupling, the high energy limits of the Green's function and of the collinear 
splitting functions are regulated by different indices, which both originate from the frozen 
coupling hard Pomeron exponent. We shall define the index uJsit) by (cfr. Sec. 14.41) 

G{Y; t, to) ~ = exp[Ljg( ~\.'^ )Y + diffusion corrections] , = ctg — - (6) 

\/2irasX"y 2 TT 

in the limit uJs{t)Y S> 1 and t ~ S> 1, and the index ujcii) by 

xP{as{k^),x) x-'^=Wp(as) , (7) 

where P{as{k^),x) is the resummed gluon-gluon splitting function (Sec. El). The exponent 
LOg in Eq. (jH)) used to be defined as the location of the anomalous dimension singularity in the 
saddle point approximation. It is now understood see also [TSj, that this singularity is 
actually an artefact of the saddle point approximation, and that the true anomalous dimension 
singularity, located at u; = uJc{t), causes the power behavior of the effective splitting function. 
This result has then been confirmed in the alternative resummation procedures of ^36137111^. 

Even the definition in Eq. © is not free of ambiguities, due to the occurrence of diffusion 
corrections to the exponent j24ll25ll51l26j which rapidly increase with Y, and to the contamina- 
tion of the non-perturbative Pomeron, which dominates above some critical rapidity [221123] • 

In the following, both regimes t ^ to and t ^ to will be discussed in detail in the RG- 
improved approach, by emphasizing our perturbative predictions and their range of validity. 

2.2 Scale changing transformations 

Let us note that the symmetrical scale choice z/q = ^^o performed in Eq. ^ is not the only 
possible one, and is physically justified only in the case k ^ ko- This configuration occurs for 
example in the process of 7*7* scattering at high energy with comparable virtualities of both 
photons [20], forward jet/vr'^ production in DIS jSj, or production of 2 hard jets at hadron 
colliders ^H]- However, in the typical deep inelastic situation, when one of the scales is much 
larger, k ^ ko {ko ^ k) the correct Bjorken variable is rather k'^/s (k^/s). In order to switch 
to this asymmetric case one should perform a similarity transformation on the gluon Green's 
function of the form 
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where A:> = max(A;, fco) ^-nd A;< = min(A;,A;o). The transformation (jSJ imphes the fohowing 
change of kernel /C^^ 

ICUk,k') ^ IC::,{k,k') = ICUk,k')(^^^ , uo = k\ (9a) 

/C^(A:,A:') ^ /CL(A:,A:') = /C^(^,A:')f-j , = k'\ (9b) 

where now /C" (/C[;) means the kernel for the upper-A:^ (lower-A;'^^ energy scale choice. 

Our goal is to find a resummed prescription for lC^{k,k') which takes into account the 
large Y terms and is consistent with renormalisation group equations. The kernel lC^{k^k') 
is not scale invariant, and it can be expanded in powers of the coupling constant as follows 

oo 

^^{k, k') = ^[a,(A:2)]"+i /C^^ (A;, k') . (10) 

n=0 

where 

= 61og(WA^) ' '=T^-^' ^''^ 

and the coefficient kernels lQ^{k^k') are now scale invariant, and additionally carry some 
cj-dependence. We shall now see how the renormalisation group constraints on /C^ and K}^ 
determine the collinear behavior of /Cj^ . 

2.3 Renormalisation group constraints and shift of 7 poles 

It is important to notice that the w-dependence of the scale invariant kernels fC^, present in 
Eq. (|1U() . is not negligible (even for the small u) values being considered) and follows from the 
requirement that collinear singularities have to be single logarithmic in both regimes k ^ kg 
and ko ^ k. If A; » fco, it is simplest to discuss the kernel in its form IC"!^, Eq. (jU^). A 
leading-log k'^ analysis for k ^ k' shows that its collinear singularities are determined by the 
non-singular part (in lo space), Ai{u>), of the gluon anomalous dimension, 

asAiiu;)=jggiio)-—, (12) 

and 

^iH = -§ + OM, {Nf = 0), (13) 

In contrast the singular part as/uj is accounted for by the iteration of the BFKL equation 
itself. 

To be precise, one has 

}CZ{k,k') ^fipexp jJ'*d(logK2) A^{uj)as{K') = (l - bas{k')log^^ " , 

(14) 

where t = log Aj^/Aq^^,, indeed showing single logarithmic scaling violations. A similar 
reasoning, yields the collinear behavior of /Cj^ from Eq. ^jp) with the opposite strong ordering 
behavior k' ^ k, which is relevant in the regime kf) ^ k. 
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But /C^ and /C^ are related to K,^ by the w-dependent similarity transformations (int-Eh) 
so that the latter must have the following collinear structure 



1 



A;2 V A; 



k' 



o-sik"^ 



as 



W2\ 



b 



1 fk 

k^ If 



e{k - k') + 

as{k^) 



a. 



b 



e{k' - k) 



■ (15) 



In this expression one can see that the w-dependence provided \ T^j '^^ essential, because 
ky/k^ can be a large parameter. We also keep the w-dependence in Ai{uj)^ in order to take 
into account the full one-loop anomalous dimension. 

By expanding in hots the renormalisation group logarithms present in the collinear behav- 
ior of Eqs. (|14|15|) . we obtain the leading collinear singularities of the coefficient kernels tC^ 
in Eq. H1U|) . This implies that, in 7-space, the corresponding eigenvalues have the following 
structure 



l-Ai{Ai + b)---[Ai + {n 



(7 + 1)"+^ 



h)Ai • • • [^1 - nb] 



1 



(16) 



1) 



where the uj dependence of Ai is left implicit. Therefore the position of the 7 ^ (7 
poles is shifted by — ^ (+^) for the kernel (|15)) with symmetrical scale choice vq = 
Through this shift one is able to resum |^ the higher order 7-poles of the kernel that are due 
to scale changing effects. 

In fact, the leading and next-to-leading eigenvalues corresponding to this symmetrical 
choice of scale have the collinear behavior 



xr(7) ^ 



1 



7+2 



+ 



1 



(7 + f)^ 



(1 



1-7 + f 
Ai{oj)-b 



(17) 



Now, in order to obtain the NLL coefficient jJJ in the expansion one has to expand 
in CO the term Xo (7) to first order with subsequent identification to a.,yn^°, and add the 
Xi 



,uj ^g]-]2is. The result for the NLL eigenvalue in the collinear approximation then reads 



(7) 



UJ 

— UJ / \ AO 1 UJ 



1.^=0 



1 

'2^ 



1 



2(1-7): 



, ^i(O) , ^i(O) 



T 



(1-7)^ 



+ 



(18) 



We note that the a;-dependent shift has generated cubic poles -.3 which seem to imply 



7^ ' (1-7)^ 

double logs log -j4-, but are actually needed with the choice of scale kk^ in order to recover 



the correct Bjorken variable k'^/s. The collinear terms with Ai{lu) have instead generated 



which correspond to single logs, log^. 



double poles ^, -^3;^ 

The double and cubic poles at 7 = and 7 = 1 so obtained are precisely those of the full 
NLL BFKL kernel eigenvalue. In fact Eq. (|18jl is a collinear approximation to the full NLL 
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BFKL kernel eigenvalue which has the following form 



(l + 27)(3-27) 



+ Ut: - TTT Xo(7) + ;tC(3) + 



36 12 y 2^ ' 4sin7r7 

f;(-ir f^^"^l^!L"^^'^ + ^^%:: . . (i9) 



n=0 



(n + 7)2 (n + 1 — 7)2 



It turns out that the collinear approximation (|18|) above reproduces the exact eigenvalue 
((T^ up to 7% [m05] accuracy when 7 g]0, 1[. This suggests that the collinear terms are the 
dominant contributions in the NLL kernel. 

In the following, we shall normally incorporate the shift of 7-poles in the form 

X^(7) = X^l(7 + f ) + - 7 + ^) , (20) 

where XnL (Xnij) have only 7^— ^(7^! + ^) singularities of the type in Eq. (fTO]) . In this 
way the collinear singularities are single logarithmic in both limits k ^ ko and ko ^ k, and 
the energy scale dependent terms are automatically resummed. The modified leading-order 
eigenvalue that we adopt has the following structure (compare (|T7jl ): 

= 2V'(1) - ^(7 + f ) - V(l - 7 + f ) , (21) 

in the case of symmetric choice of energy scale vq = kko. This form of the kernel was 
considered previously in j39U4Uj . It is obtained from the leading order BFKL kernel by 
imposing the so-called kinematical (or consistency) constraint |4 H I42 | B3] which limits the 
virtualities of the transverse momenta of the gluons in the real emission part of the kernel. The 
origin of this constraint is the requirement that in the multi-Regge kinematics the virtualities 
of the exchanged gluons be dominated by their transverse parts. The NLL contribution of 
the resummed kernel, Xi was then constructed by the requirement that the collinear 
limit in Eq. (|17|) should be correctly reproduced, and the exact form of the NL kernel ()19|1 
should be obtained also. 

The final NLL eigenvalue function proposed in |lUllllj reads 

1 

Xiil) = Xi (7) + 7^x0(7)— ^ — 

2 sm TTj 
-A^m'{J)-[A,{0)-b]^l^'{l-J) 
+Ai(a;)^'(7 + f ) + [Ai(a;) - 6]^'(1 - 7 + f ) 

-Y[xo(7)-X^(7)]. (22) 

The first line is the original NLL term Xiil) with the subtraction of the cubic poles which 
come from the changes of the energy scale and which are resummed by the leading order 
a;-dependent kernel ((22) • The second and third lines contain shifted collinear double poles, 
and finally the last line contains the shifted single poles which additionally appear as an 
artefact of the resummation procedure. 
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2.4 cj-expansion and collinear resummation 

In the present paper we choose a form of the improved kernel that differs somewhat from 
that of Ref. ^] — quoted in Eqs. 1)21122(1 — by using the possibihty of translating part of 
the a^-dependence in Eq. (fTU]) into additional w-dependence. Actually, it was pointed out 
in jlUUllj that, at high energies, w is a more useful expansion parameter than as{k'^), the 
relation being given roughly by u; ~ asXOy as noticed already in connection with Eq. ()18|) . 
The w-expansion is a systematic way of solving the homogeneous equation 

{io - JC^)J'Uk) = , (23) 

where /C^ is given by Eq. (jTU)) . by the 7-representation 

^.(t) = /^e^*-i^^"W, (24) 

in which Xioil) = ^Lil) satisfies a non-linear integro-differential equation equivalent to 
Eq. (|23j) . The latter is derived by using the representation t ^ —d-y in Eq. ()lfl|) . and is given 

by d 

\ / \ 2 



Approximate solutions to Eq. ()25() can be obtained either by truncating at, say, NL level (i.e., 
setting X2 = X3 = ■ ■ ■ = 0); or by expanding in the w-parameter to all orders. The latter 
procedure yields the solution [TUIITT] 



= Xo (7) + ^—J77X + ^ 



X^(7) X^(7) 



xni)^^fxf\' fx^ 



x^(7) vx^y VX^ 



+ ..., (26) 



and amounts to replacing the kernel )Cuj by an effective kernel as(/c^)/CJjf , where /C°^ is scale 
invariant. The corresponding characteristic function Xujil) Eq. (|26() is — very roughly — 
obtained by the replacement — > tj/xo; in Eq. (|Tn|) . so that indeed uj plays the role of a 
new expansion parameter. A virtue of the expansion (|26|) is that it contains simple (leading) 
collinear poles only, because the double-poles left in Xi after the w-shift are canceled by the 
denominators. 

The w-expansion is particularly useful for the resummation of the leading collinear sin- 
gularities of Eqs. (fT3)) and (|16() . Suppose we first take ctg frozen (limit 6 = 0). Then, the 
leading poles of Eq. (|16|) have approximately the factorized form 

X--X;^(X-)", X- = ^ + j^^, (27) 

7+2 1-7+2 

(valid for 7 + ^~0orl — 7 + ^ — 0), so that the resummed behavior (|15() reads 

00 

ICu^^Yl (^sK^[asKr = asK^il - asK'^)~^ . (28) 

n=0 

Exactly the same result can be obtained by the w-expansion H26|) truncated at the NL level, 
by setting Xi/Xq — Xc > ^.nd thus considering the kernel 

iC^ = as{K'^ + u:K^) . (29) 



9 



In fact, the resolvent of the latter is given by 

a. ^ - K:^]-^ ^{l- a.K^r^ [u: - asK^{l - a^K)-^] , (30) 

and is then proportional to the Green's function of the resummed kernel (|28j) . 

In other words, leading-log collinear singularities are equivalently incorporated by a string 
of subleading kernels (as in Eq. (|28|) ). or by a NL contribution of order cksUJ (as in Eq. (|29]) ) — 
apart from a redefinition of the impact factors. In the realistic case with running coupling it is 
straightforward to check that 6-dependence only remains in the first term of the w-expansion 
(EHl) 

^ + a; (-4^ + . ^7^^ ) + ■ • • , (31) 
V7+2 1-7+2/ 

whereas it cancels out in all remaining subleading terms. Therefore, in order to incorporate 
the leading log collinear behavior in the form (|31|) we can set, for instance, 

= as{q^)K'^ + ujas{kl)K'^ + NLL , (32) 

as an improved leading kernel. Here we assume that the scale for otg in the leading BFKL part 
is provided by the momentum of the emitted gluon q = k—k', as suggested by the 6-dependent 
part of the NLL eigenvalue in Eq. ((T^. which corresponds to the kernel b-ij log ^ | j^^g (see |2]), 
and — via cj-expansion — to the 6-term in Eq. pif) . A simplified version of Eq. ()32() without 
the NLL term and with one collinear term (for 7^0) was used in for a phenomenological 
analysis of the structure functions. 

Note that, if we take literally the w-expansion (|26jl with the choice of NLL term (|22|) . 
then Xi/Xo would coincide with Xc close to the collinear poles, but would be different in 
detail away from them, and would actually contain spurious poles at complex values of 7 due 
to the zeroes of Xoil)- Such poles cancel out if the full w-expansion series (|2H|) is summed 
up, but are present at any finite truncation of the series, thus implying poor convergence of 
the solution whenever 7-values close to the spurious poles become important. For this reason 
in this paper we prefer to resum collinear singularities by the improved kernel ()32() . which 
contains only collinear poles. Furthermore, the NLL term needed to complete Eq. (|32|1 — 
to be detailed in the next section — turns out to have only simple (leading) collinear poles, 
because the running coupling terms have been already included in the g^-scale dependence 
of the running coupling. Therefore, the full kernel has the same virtues as Eq. (|26p in the 
collinear limit and, lacking spurious poles, is more suitable for numerical iteration. 

3 Form of the resummed kernel 
3.1 Next-to-leading coefRcient kernel 

We have still to incorporate in our improved kernel the exact form of the NLL result 
in the scheme of the ots expansion, i.e. 1)32(1 . We choose to start from the leading kernel in 
Eq. I|32|) which incorporates both the collinear resummation and the running coupling effects 
due to the choice of scale q^. The full improved kernel then has the form 

IC^ = as{q^)K^ + iJas{kl)K^ + al{kl)kt , (33) 

where A;> = max (A;, k'), fc< = min(A;, k'), and K'^ is determined below. 
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We recall that the Mellin transform of the colhnear part /C^, defined by 



leads to the expression 



/C(*,t') = 4Pf|i) . (35) 



One can match the above prescription to the standard kernel at NLL order by expanding in 
UJ and in bois to first order 

K:^ ~ a,(fc2)(i^o + coK^ + ujK^) + a2(^o + i^--), (36) 
where we have defined 

b 



K^ = KT\ K^ = Kr\ K}- ' 



xri7) = -^^{Xo + xl), (37) 



by noting that the running coupling term has the form [see Eqs. (|88I89() and Add. \K\ 



(38) 



Reg 



By replacing the expression (|36jl into Eq. we obtain the relationship with the customary 
BFKL Green's function 



UJ 



IC^Y^ = (1 - as{Kl + i^°)) ^[u-as (Kq + a^i^i + ©(a^))] ^ , (39) 



where Ko and Ki are LL and NLL ^-independent kernels. The two expressions will match 
provided we identify 

Ko = Kl 

kl = K,- K^,{Kl + K") - Kr , (40) 

and we properly redefine the (so far unspecified) impact factors (see Sec. 1^1). Thus the term 
in (|4U() corresponds to the customary NLL expression (|19() with subtractions. 
Li 7-space the subtracted NLL eigenvalue function which corresponds to the Kf has the 

following form: 

Xi(7) = Xi(7) - X'omxhh) + X°(7)] - Xr(7) 

= Xi(7) + 2>^o(7)-^ - ^o(7):^(^ + 2 (X^ + xg) • (41) 

The subtractions cancel the triple poles (due to change of energy scales) and the double poles 
(from the non-singular part of the anomalous dimension). Therefore the resulting kernel xi 
contains at most single poles at 7 = 0, 1. Eq. together with the eigenvalues (PT|). (|5Hl 
and ()41() gives a complete prescription for the resummed model. This new formulation is 
identical to the previous cj-expansion |l(JUllj near the collinear poles. It has the advantage 
that it can be easily transformed into the (x, fc^) space (it is free of ratios in 7-space, such as 
Xi/Xo) aiid avoids the spurious poles that were present in (|26l) . 
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Note that the choice of scale in a<j in the first term in Eq. H33|) is determined by the form 
of the NLL part. Any change of scale in this term would correspond to the change of NLL 
terms proportional to b. The scale for the coUinear parts is chosen to match the standard 
DGLAP formulation whereas in the NLL part is purely conventional, and its change would 
be of the NNLL order. In the following, in order to study the dependence on renormalisation 
scale uncertainties, we introduce the quantity and generalize eq. ()33() as follows 

IC^ = {asixW) + hal log xl) + u: (a,(x2 A;^ ) + ha^ log xl) + ^^(x^ fc^ . (42) 



3.2 Form of the kernel in {x, k^) space 

We define the resummed kernel in {x,k'^) space as the (integrated) inverse Mellin transform 
of JCoj-. 

iC{z; k, k')= [ ^ z-^ -iC^{k, k') (43) 



27ri uj 

where the real variable z can assume values between x and 1. 

The subtractions of (|lT|l are translated into (x,A;^) space to give 



\xo{l)—^, — r 
I sm (7r7 j 



1 



log2^ + 4Li2(l-^ 



-Xo(7) 



^i(O) 
7(1-7) 



4|fc2 - A;'2 
-^i(O) sign(A:2 



1 , \k^ 



Ml 



Q 



1 1 1 



1 , \k'^ 



fc2 



(44) 



where the dilogarithm function is defined to be 

dt 



Li2(tt;) := 







t 



log(l-i), Li2(l) 



TT 



(45) 



In (x, /c^) space the symmetric shift is translated into the symmetric kinematical constraint 
which has to be imposed onto the real emission part of the BFKL and also into the collinear 
non-singular DGLAP terms: 

(46) 



k 

kz < k' < - 

z 



(in the following we denote the imposition of the kinematical constraint onto the appropriate 
parts of the kernel by the superscript (kc), i.e. KQ^{k, k')). 

The final resummed kernel /C(z; k, k') is the sum of three contributions: 



1 dz 

z 



dk'^ iC{z;k,k')f{-,k') 

z 



1 dz 
z 



dk 



12 



as{q')K{z- k, k') + as(ki)K^'{z- k, k') + ai{ki)Ki{k, k')^ /(-, k') . 



(47) 



12 



The different terms are as follows: 

• LO BFKL with running coupling and consistency constraint (q = k — k') 



^ Idk 



12 



as{q^)K^'{z;k,k')\ f{-,k') 



dz f d^q 
irq'^ 



/(-, \k + q\)Q{- - k')Q{k' - kz) - Q{k - q)f{-, k) 



non-singular DGLAP terms with consistency constraint 

"1 dz 



dk" a,{ki)K^^{z-k,k')f{-,k') 



1 dz 

z 



dk'' 



k ^ 



{kzY 



as{k')z-P,g{z-)f{-,k') 



+ — 



Z fc'^ 

NLL part of the BFKL with subtractions included 



dz r^^l^) dk'^ _ 2. k' ~ k' X 

as{k )z-Pgg{zj)f{-,k) 



dz 



z 



- / dk'^ al{kl)K^{k,k')f{-,k') 



z 



47. ~J 



dk'^ al{kl) 



V 9 



3 J |A;'2-A;2| 



+ 



1 



2 2 

+ TIT + 



32 V fc2 \k'^ A;2 



1 1 



log 



A;2\\ AU2{l-kl/k 
pjj + \k'^-k^\ 

1 



-4^i(0)sgn(A;2 _ k'^) 



2 



\k'^-k'' 



1 |fc^^-fc2| . 



3.^ 



fc^' + fc2)2 

32fc'2A;2 



1 



vr 



fc2 + y2/j/2 



log 

1 - y 



^ J X Z Z 



The non-singular splitting function in the DGLAP terms is defined as follows: 

P -P -i 

^99 ~ ^99 , ' 



(48) 



(49) 



(50) 



(51) 



where we take 



I- z 



99 



+ z{l - z) + 



Z 11 c-/-, N 

^ \ (5(1 - z) , 

(1 - z)+ 12 ^ ^ ' 



(52) 



(we only consider purely gluonic channel, nj = 0). Also we note that the argument of the 
splitting function P has to be shifted in (|49|1 in order to reproduce the correct collinear limit 
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when the kinematic constraint {kz < fc' < |) is included. This follows from the inverse Mellin 
transform of Eq. H35|) 



K^^{z;k,k') 



duj Ai{(jj) 
kl 



A;. 



kl 



P 



(53) 



In other words the correct variable in the splitting function is modified by the ratio of two 
virtualities in the case when the kinematical constraint is included 



k 

z— < 1 for k' < k 
k' 

k' 

z— <1 ioi k < k' . 
k 



(54) 



3.3 Choice of scheme 



The prescription formulated above for the kernel eigenvalue (|4ip is free of double and cubic 
poles in 7 = (and 7 = 1), however there are still some residual single poles. These poles 
come from the constant terms from the expansion of subtraction Xo + Xc around 7 = 
(7 = 1). Expanding this subtraction around 7 = one obtains 



-xl 



1 

2^ 



7 



therefore there appear additional singular terms, 

^2 



TT 



^i(O) 



1 

7 



(55) 



(56) 



in the subtracted kernel x? which are not shifted. Furthermore, the term (|56jl contributes to 
the 2-loop anomalous dimension, together with the constant term arising from the leading 
kernel as follows: 

1 +u)Ai 



where 



Xo + ^Xc 

C{u) 
C(0) 



7 + 



a.C(a;) + 0(7 + f), 



(57) 



^i(c^) , V(i + w)-V'(i) 



+ 



Am . 

By combining (|56() with 1)571) we would get the contribution 

sCiuh^^^ ~ ^ [C(0) - Ciu) (1 + coAiiio))] 



A7(2) = ^c(O) 



(58) 



(59) 



where s(l + ujAi(iv))/uj is the DGLAP anomalous dimension in the leading order. 

The expression ()59() violates the momentum sum rule A7(2)(u; = 1) = 0. 

We thus consider two possible forms of subtraction. In the first scheme A we add and 
subtract from the NLL part the term proportional to C(0) in the following way. 



Xi(7) - xr(7) = Xi(7) - C(0)xo(7) + C{0)x^{^) , 



(60) 
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which leads to the following modification of the kernel in (x, A;^) space 



— / dk 



' '"^{as{q^)K^'{z;k,k')+as{kl)K^^{z;k,k') + 

(kl) [Ki{k,k') + C{0)K^'iz;k,k')-CiO)Ko{k,k')\}f{^,k'). (61) 



This scheme satisfies general RG constraints, but contains the anomalous dimension H59|) and 
violates the momentum sum rule. 

In the second scheme B we shall consider a modification which adds the shifted pole to 
the NLL kernel with the oJ-dependent coefficient (1 + u;Ai(uj)) 



1 1 

- + 



7 1-7 



1 



1 



\7"r2 -'■"'"2 7/ 

(62) 

It is straightforward to check that in this case the 2-Ioop anomalous dimension vanishes^, 
due to a cancellation between the pole term (|62|) and the constant term in ()57() . Therefore, 
scheme B satisfies energy-momentum conservation. 

The change in the resummed kernel in (x, A;^) space corresponding to scheme B is obtained 
by inverse Mellin transform of ()62() and is given by 



1 dz 

z 



dk'^ {a,(q2)i^kc(^. ^ asikDK'^^iz- k, k') + al{kl)ki{k, k')] /(^, k')- 



al{e)fC-,k') + 



°° dk'"^ 



,x 



^2 K'^ Z 



k^ 



k' 



fc2 



-^al{k'')zjS,{z'-)fC-^k') 



k\ 



(63) 



with the function 6*2(2;) given by 



S2{Z) 



lUz 



132 + 24tt^ + z[-541 + 247r2 + 72z(l + 3z)] - 1441og(-l + i) log(^) 



+ 12{ log(l - z)[-l - 2z(23 + zi-15 + 8z)) - 12(1 + z) log(l - z)] + 12z log(-l + i) log(^) 

+ 2z[l + z(-21 + 5z)-6 log(l - z)] log(z) - 6(-l + 2z) log'^{z) 

2 



+ 144(-l + z) [Li2(z) + ^log(i)log 



:i-2)2 



vr 



y]- 144(1 +2z)Li2(l-z)!> . (64) 



Note that whatever scheme we choose, Ki contains higher-twist poles (at 7 = — 1, —2, . . . 
and 7 = 2,3, . . .), which are not shifted. In the calculations that follow we keep these poles 
unshifted independently of the choice of energy-scale. This means that calculations of the 
Green's function carried out with different energy-scale choices will formally differ at NNLL 
level. In practice however we find that this energy-scale dependence is very small. 

^We use here a generalization of the Qo-scheme |44l . We do not try to include the known 2-loop expression 
in the MS scheme because it is subject to a scheme change and to kernel ambiguities which are not fully 
understood yet. 
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4 Characteristic features of the resummed Green's function 



We shall first investigate the features of the two-scale Green's function^ G{Y; k'^, k^) based on 
the form of the resummed kernel just proposed. In the perturbative regime fc^, /cq ^qqd 
with iOs{k'^)Y large we have both perturbative contributions, leading to the hard Pomeron 
exponent, and non-perturbative ones, due to the asymptotic Pomeron, which is sensitive 
to the strong coupling region. It was noticed in |2HI22j that the hard Pomeron dominates 
for energies below a certain threshold as{k'^)Y < l/bcop beyond which there is a tunneling 
transition to the non-perturbative regime. It has also been noticed j22j, that in the formal 
limit 6 — > with as{k'^) fixed the Pomeron is suppressed as exp(— l/^a^), so that one can 
define a purely perturbative Green's functions and investigate the diffusion corrections to the 
hard Pomeron exponent. In the following, we use the 6-expansion up to second order, so as to 
obtain the exponent uJs{t) and the additional parameters occurring in the diffusion corrections 
predicted by our improved small- j; equation. Furthermore, we analyze the perturbative non- 
perturbative interface numerically so as to estimate, as a function of logQ^, the critical 
rapidity beyond which the non-perturbative Pomeron takes over. 

Since the perturbative rapidity range turns out to be considerably extended with respect 
to LL expectations, we shall be able to extract numerically the full perturbative Green's 
function and among other things its high-energy exponent and diffusion corrections to it. 



4.1 Frozen coupling features 

Let us first consider the features of G{Y;ti,t2) in the limit of frozen coupling = as(^o)' 
i.e. 6 = 0. In such a case the kernel /C^ becomes scale invariant, but the solution to Eq. Q 
is still non-trivial, due to the w-dependence which complicates the y-evolution, it no longer 
being purely diffusive. In fact, the characteristic function becomes 

asXw{7,as) = Usixo +uJXc) +alxi , (65) 

and the important to values, corresponding to the pole of the resolvent, are defined by 

w = asXuj{'y,as) , (66) 

whose solution at fixed 7 we denote by 

w = "sxiff (7>as) , (67) 

the superscript (0) referring to the 6 = limit. The effective characteristic function (|67)1 
so defined has the interpretation of a BFKL-type eigenvalue reproducing the pole (|^^ . As 
such, it can be compared, at least for frozen coupling, to the analogous quantity defined in the 
"duality" approach of Ref. pEj. It provides information about the hard Pomeron exponent 
and the diffusion coefficient D = Xm/'^Xm- In Fig. ^we compare the results for the exponent 
ujs as a function of calculated in the case of fixed coupling for schemes A, B and the original 
w-expansion method presented in jlOllllj . The critical exponent is obtained by evaluating 
the effective kernel eigenvalue at the minimum 

^ = asxS (7-, «s) . (68) 

^In Sees. |1| and |K| we remove for simplicity the ~ symbols used before to denote RGI quantities in our 
present scheme. 
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Figure 1: ujg as a function of as for different subtraction schemes together with the original 
result for the uj-expansion. The calculation is done in the fixed coupling case. 



All resummed results for the intercept are significantly reduced in comparison with the LL 
result and they all give stable predictions even for large values of otg- As we see from the 
plot, the changes of resummation procedure as well as subtraction scheme do not significantly 
influence the values of iVg- They give at most 20% change at the highest a<j ~ 0.35. In Fig. |21 
we show the effective kernel eigenvalue as a function of 7. We have considered here the 
asymmetric w-shift, which corresponds to the upper energy scale choice z/q = k"^. In this case 
it is easy to show that close to 7 = the effective eigenvalues from scheme B and the original 
(j-expansion [TJ satisfy the momentum sum rule. This is illustrated in Fig.[2lby the fact that 
fisXeff(7 = 0,q;s) = 1 for all values of Oig in these schemes. This can be seen by expanding 
around 7 = 0, where we have 

Xuj[l,as)cc (69) 

7 

which for 7 = gives ojAi[ijo) = — 1, which has the solution uj = 1. Note that a second fixed 
intersection point of curves with different as occurs at 7 = 2. This is expected from energy- 
momentum conservation^ in the collinear regime Qq ^ Q^, because of a behavior similar to 
Eq. H69() around the shifted pole 1 + a; — 7 = 0. This intersection has no counterpart in the 
approach of Ref . • 

We also examine the second derivative Xesil^^s) which controls the diffusion properties 
of the small-x equation, Fig. |21 As we see from the plot, the second derivative is more model- 
dependent than the intercept tOs, though the two models A and B presented in this paper 

■^Such an intersection occurs in scheme A also (where momentum conservation is not satisfied) as an artefact 
of the collision of the shifted pole at 7 = 1 + a; with the unshifted one at 7 = 2. 
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Figure 2: asXefr(7, Os) as a function of 'j in different schemes for different values of ag: 
as = 0.1 (dash-dotted line), Ug = 0.2 (solid line), Ug = 0.3 (dashed line). The calculation is 
done in the fixed coupling case. 



give quite similar answers. The value of the second derivative will influence the diffusion 
corrections to the hard Pomeron, as we shall see in Sec. 14.41 and also the transition of the 
solution to the non-perturbative regime. 

4.2 Numerical methods for solution 

In this section we are going to investigate in detail the shape of the solutions to the integral 
equation^ with the resummed kernel given in sections 13.21 and 13.31 To this aim we solve 
numerically the following integral equation^ 

G{Y- k,ko) = {k, ko)e{Y)+ dy dk'^ IC(Y - y; k, k')G{y; k' , ko) (70) 

*An interesting iterative method of solution to the NLL BFKL equation has been recently proposed (45| . 
By using this method it is possible to solve the equation directly in [x, k) space and keep the full angular 
dependence. 

®Here we change slightly the notation in the first argument of /C, writing log ^ — Y — y instead of z. 
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Figure 3: x'^lmjCts) as a function of as for two different subtraction models and the uo- 
expansion scheme. 

with (so as to have the same normahzation as in Eq. ©), 

27rfe2G(o)(fc,fco) = <5(log^). 

We use the method of iterations and discretized kernel similar to that introduced in j33j . 
More precisely in our problem (see (|47j) ) we can rewrite the kernel in the following way: 

lC{Y-y- k, k') = J2 ^a{Y-y- k, k') = J] K^ik, k') P^iY-y) Qr (y-y-max(log p log |)) , 

a a 

(71) 

where the index a enumerates different terms in the equation H47() (that is LL BFKL, LL 
DGLAP, and the different components of NLL BFKL with subtractions), each of which 
factorize into transverse and longitudinal parts. The Pa are the singular and non-singular 
pieces of the splitting function as well as the subtraction terms 52 (x). The additional 0r 
stands for the kinematical constraint, applied to all terms that in Mellin-space have an w-shift. 

To find the solution numerically one introduces a grid in rapidity Y and logarithm of 
momentum, r = logk/kQ, with small spacings. Ay and At respectively. The solution is 
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then calculated at the grid points. Linear interpolation gives the values of the solution in the 
points between the nodes of the grid 

G{Y; k, ko) = Y,Yl G{Yi; kj,ko) (72) 

i j 

where (j)i{Y) and ipj{k) are the appropriate basis functions for linear interpolation. To find 
the solution for G the equation (jTUf) is solved by a method of evolution in rapidity. In a first 
step one takes G{Yq = 0; km, ^o) = G^^\km, ko) and estimates G{Yi; km, ^o) at the next point 
of the grid, Yi = AY, using the integral equation (|7flj) . This gives a first approximated value 
for G{Yi;km,ko). This function is then again used in equation (jTUj) to calculate the next 
approximation. Usually a few iterations are sufficient to find an accurate answer (typically 
5 — 8). After obtaining G(Yi; km, ko) with the desired accuracy one proceeds to calculate the 
solution on the next point of the grid Y2 = 2AY and so on. The procedure is then repeated 
for all points of the grid in rapidity Yn = nAY. 

The procedure presented above requires numerous evaluations of the right hand side of 
equation (|7()j) . Given the fact that we have two convolutions in y and k' , with the complicated 
kernel /C, such a procedure can be quite time consuming. 

In order to speed up the calculation one can discretize in k' the kernels and in y the 
functions Pa using the basis functions in the following form 

i^ii = j dk'^Ka{km,k')Mk') 

Pi"Jj = j dyPa{Yn-y)<^j{y) , (73) 

where we have used the fact that the functions Pa depend only on the difference Yn — y which 
— together with the linear interpolation — results in a one-dimensional vector P instead of 
a matrix. One can simplify the treatment of the Qr function in (|71j) by using the same grid 
spacing in y and in log A;, Ay = At (or for energy-scale choice fo = k'^. Ay = 2 At). After 
the discretization procedure, the convolution on the right hand side in equation (|7flj) (and 
using (|71|)) can be then represented as a multiplication as follows 

„ „ iniax n-\m-i\ 

dy dk'^ ICaiYn-y; km, k') G{y- k' , ko) = Y, Yl Pn-j-\m-^\ GiVv ko) , (74) 

i=0 j=0 

so that in practice all the integrations present in equations ()7U() are performed once before the 
evolution, and then only the multiplications of kernel matrices and gluon Green's function 
vectors are done during the iterations. 

Of course, in a numerical analysis one is not able to use exact distributions — in particular 
for the delta function in k as an initial condition, see Eq. Q- In practice what is done is to 
set to 1/Ar one point on the fine grid i.e. 27rkQG^^\km, ko) = -^Smo, where At is the grid 
spacing in log k. The resulting Green's function will be finite in the Y = limit but dependent 
on the size of the grid spacing. We illustrate this effect in the upper set of curves of Fig. |^ 
where we have solved the equation (|7()|) with the kernel in LL approximation with 3 different 
grid spacings At = 0.05, 0.1, 0.2. One might be worried by the apparently substantial 
dependence on the choice of the grid spacing At. However this is just a consequence of the 
grid-dependent discretization of the initial 5-function and disappears when one convolutes 
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Figure 4: Gluon Green's function as a function of rapidity Y for three different grid spacings 
At = 0.05,0.1,0.2. LL evolution is used with a fixed coupling, ag = 0.2; e ~ 0.2/co. 

the gluon Green's function with some smooth impact factor. We will therefore consider from 
now on a slightly asymmetric choice of scales, G{Y;to,to + e) with e = 0.2. In the lower set 
of curves of Fig. |1] one sees that the dependence on the grid spacing in this case is relatively 
small. For the remaining plots in this paper we have used At = 0.1 or smaller. 

4.3 Basic features of the Green's function 

Let us now discuss the properties of the gluon Green's function obtained with the method 
discussed above. We shall use a one-loop coupling with n/ = 4, normalized such that 
as(9 GeV^) = 0.244.^ The coupling is cut off at scale k = 0.74 GeV — a detailed anal- 
ysis of the sensitivity to this regularization is postponed to section 14.51 In the kernel, for 
the time being we consider Uf = 0, since our single-channel RGI approach does not properly 
account for the quark sector (however we will see below that simply varying uj in the kernel 
has only a small effect). 

Results will be given given for: LL evolution (with as(g^)); our two resummation schemes, 
A and B; and two variants of 'pure' NLL evolution: one, labeled 'NLL as{q'^)^ where the kernel 
is as{q'^)KQ + a^{k^)K^^^ , with K^^^ corresponding to eq. (fTU]) without the first term in 
square brackets; and another, labeled 'NLL as{k'^)\ where the kernel is asik'^)Ko + a^{k'^)Ki, 
and Ki corresponds to eq. (fT^ in full. 

Fig.Elshows Green's functions G{Y; ko + e, ko) as a function of rapidity Y, and fig.Elshows 
kkoG(Y; k, ko) as a function of k ior Y = 10. To aid legibility, each figure has been separated 
into two plots, the left-hand one (a) showing LL and schemes A and B, while the right-hand 
one (b) shows the two pure NLL curves and scheme B. We choose a moderately high value 
for the initial transverse scale, ko = 20 GeV, as{ko) ~ 0.15, so as to be able to focus on 

®As one obtains, roughly, by running as{M^) = 0.118 down to 9 GeV^, taking into account flavor thresholds 
and the two-loop /3-function. 
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Figure 5: Gluon Green's function G{Y; ko + e, ko) as a function of rapidity Y : (a) for LL and 
the two RGI schemes A and B; (h) for scheme B and two variants of pure NLL evolution. 
The parameters are ko = 20 GeV and e ~ 0.2ko. 



the perturbative aspects of the problem (non-perturbative effects are formally suppressed by 
powers of A^/feg). Such a scale has been used for BFKL dijet studies at the Tevatron 

A number of features of fig. are worth commenting. Most noticeable is the significant 
reduction in the high-energy growth of the Green's function when going from LL evolution 
to our resummed schemes A and B. This is as expected from the discussion of high-energy 
exponents, fig-E Also important is the fact that for the RGI schemes the high-energy growth 
does not start until a rapidity of about 4. Both of these observations are relevant to the 
problem of trying to reconcile theoretical predictions with the lack of experimental evidence 
for a strong high-energy growth of cross sections at today's energies. The small difference 
between the two RGI resummation schemes, A and B, is in accord with their slightly different 
ujs values (cf . fig. . 

As regards the transverse momentum dependence of the Green's function, fig. there 
are a number of further differences between the LL and RGI results. The higher overall 
normalization for LL evolution is just a consequence of a larger lOs value. But one also sees 
that the large-A; tails in k for the resummed models are substantially steeper than in the LL 
case. This can be understood by comparing the diffusion coefficients in these models: the 
RGI models are characterized by a smaller x'es and, as a consequence, they have less diffusion 
than in the LL case. As was the case for the Y dependence, the two RGI schemes give very 
similar results, here differing essentially only in the normalization. 

Some comments are due concerning the structure at low k: there, there is a component 
of the evolution that is sensitive to the larger coupling, 0^(1 GeV^) ~ 0.4. For the LL case 
the resulting stronger evolution (than at k^) over-compensates the suppression due to the 
large ratio of scales k^/k, leading to the absence of a decreasing low-A: tail. For the RGI 
schemes the difference between uJs values at 1 GeV and kQ is not sufficient to bring about 
this overcompensation for Y = 10, so there still is a decreasing tail for small k. However 
the results are sensitive to the fact that at large the difference between cOg values for the 
two schemes becomes non-negligible. This is what causes the low-fc Green's function to be 
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Figure 6: Gluon Green's function 2TTkkoG(Y;k,kQ) at rapidity Y = W as a function of the 
transverse scale k. The sets of kernels used in plots (a) and (h) are the same as in figure\^ 



almost three times larger for scheme A than scheme B. It should of course be kept in mind 
that all the properties at low k are strongly dependent on the particular choice of infrared 
regularization of the coupling. 

Let us now examine the right-hand plots of figures El and |H1 which show results with pure 
NLL evolution. We recall that the original motivation for introducing RGI resummation 
schemes was the large size of the NLL corrections, and in particular the fact that for moderate 
values of the coupling the NLL terms change the sign of xil) ^-iid its second derivative around 
7 = 1/2, with the situation being even worse in the collinear region. Nevertheless, as was 
pointed out by Ross 0, because of the change of sign of x"{\)^ the usual saddle point at 
7 = ^ is replaced by two saddle points off the real axis, at 7 = 1/2 + ivq and 1/2 — iv^, and 
it is the value of x &t these new saddle points that determines the high-energy behavior of 
the (fixed-coupling) NLL Green's function: 

(75) 

Since x{l) = X*(7*) this gives 
MoG{Y-k,ka) ^ e ^0 cos (^KM log + a,y 9[x(- + xi^o)] J • (76) 

When Xi(7) is symmetric in 7 <-> 1 — 7, as is the case if we use as{q^) in the LL term (or as 
can be achieved with the modified Mellin transform suggested in ^ and used in [Jj), then 
OisX{\ + ^0) is real, having a value of about 0.2. One therefore expects to find a high-energy 
growth of the Green's function that numerically is not so different from that with out RGI 
resummed schemes. This is precisely what is observed in fig. [Hb for the 'NLL as(g^)' result. 

On the other hand if Xi{l) is not symmetric in 7 <-> 1 — 7 then x will be complex at 
the saddle points. This is the case for the 'NLL agik'^Y kernel and the change in sign of 
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the Green's function around y = 18 can be understood as a direct consequence of a zero of 
eq. ^ when a^Y ^[xik + i'^o)] = 7r/2- 

The oscillatory behavior of eq. (|76p also becomes an issue when k ^ kQ, as is visible in 
fig-Eb- For NLL evolution with as{q'^) the change of sign intervenes only for ratios of k/ko 
that are fairly small or large from a phenomenological point of view (at least for Mueller- 
Navelet or 7*7* type processes). For evolution with asik"^) the situation is more dramatic 
because of the sum of terms in the argument of the cosine of eq. (|76() . 

So our overall conclusions regarding NLL evolution is that, while in certain instances 
it may give results that are not too different from those with RGI methods, in general it 
offers only limited predictive power, because of the strong sensitivity to the details of the 
formulation. Though here we have just discussed renormalisation scale sensitivity, we note 
that changing the energy scale vq from say kkQ to k'^ also leads to a Green's function that 
oscillates as a function of Y, since once again the characteristic function is asymmetric. 

A final point relating to figures and 101 concerns the overall normalization of the results. 
One sees that at low Y the LL and NLL results all have similar normalizations, while the RGI 
results are slightly lower. This is because the w-dependence is associated with an implicit 
NLO impact factor. This of course has to be taken into account should one wish to use the 
RGI Green's function in conjunction with any NLO impact factor calculation, as is discussed 
in detail in section |21 To close this section we present brief results on n j and renormalisation 
scale dependence for the RGI schemes. 

Our RGI approach has been constructed for a purely gluonic channel and only scheme 
B satisfies the momentum sum rule in this case. For phenomenological purposes one would 
wish to include quarks and in Fig. [7| we present the two schemes in the cases when n j = 
and Uf = 4 in the NLL kernel. As is clear from the plot, having nj ^ does not change 
the result in a significant way. We note that the full inclusion of quarks in a RG-consistent 
manner is a non-trivial operation in this framework especially if one is to construct resummed 
quark anomalous dimensions that satisfy the momentum sum rules. 

Finally, we show the dependence of the gluon Green's function on the renormalisation scale 
choice eq. (|I2|) . We have varied the scale in the range 1/2 < xl<2. The results of the 
calculation are presented in Fig. |H1 where the yellow bands correspond to the renormalisation 
scale variation for two resummation schemes. 



4.4 6-expansion of intercept and of diffusion coefficient 

In order to properly evaluate the hard Pomeron intercept ujs in the case with running coupling 
it is necessary to control the corrections with respect to the frozen coupling limit. To this 
end we shall apply the 6-expansion method presented in |23j . 

According to this method, we use the formal limit b ^ (with as(to) kept fixed) in order 
to suppress the non-perturbative Pomeron. The left-over perturbative Green's function can 
then be expanded in b in the form 

G{Y- to, to) = {Y- to, to) exp [bLoi^'>Y + 0{b^4Y^)j (l + 0{b^aiY^)) , (77) 

which shows a shift of uj^ of order ba^, as well as diffusion corrections of order {bas)'^{uJsY)'^ 
and {bas)'^{LJsY)'^ . The purpose of this subsection is to compute wi^^ [defined by Eq. ((77|) ] 
and the Y"^ terms both analytically and numerically. Further corrections to lOs of order b^a^ 
appear as subleading contributions in this expansion and are probably not really meaningful, 
given the complex y-dependence of the exponent involving the parameter ba^Y (23j . 
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Figure 7: Gluon Green's function G{Y; ko,ko + e) as a function of rapidity Y for two different 
resummation schemes A,B in the n/ = 0, 4 cases. Parameter e ~ 0.2/co 



We start by expanding the Os-dependence of the kernel around the frozen-couphng hmit 
up to 0(6^) by setting, for instance at scale q'^, 

(2 \ / 2 2 

log I2 + - *o)j + b^al nog2 I2 + 2(t - to) log + " *o)^ 

(78) 

where oq = cks(^o) throughout this section. We then define the kernel with frozen-coupling 
/cL°^ = JCnjl^ and the correction kernel A as 

A{t, t') ^ /C^ - = /C, - ao (i^o" + ^K) - «o^r (79) 
= Ao(t - t') + {to - t)Ai{t - t') + {to - tfA2{t - t') , (80) 

where the Aj's are scale- invariant, and are obtained from the definition (|47|) by picking up 
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Figure 8: Gluon Green's function G{Y; ko, ko + e) as a function of rapidity Y for two different 
resummation models A,B. The bands represent the uncertainty due to a variation of the 
renormalisation scale in the range 1/2 < < 2. Parameter e ~ 0.2ko 



the relevant terms in the running couphng expansions of type H78|) . We obtain: 
Ao 



log ^K^ + log ^{u^K + 2«o^r 



+ Oib'al) 



ban 



+ ujK^ + 2ao^i^ + 0{b'^ap 

d 



bal^Ki'>\ao;k\kl) + 0{b^4) 



A, 



(81a) 

(81b) 
(81c) 



Now we evaluate the Green's function G^{t,tQ) in cj-space up to second order in 6, with the 
purpose of deriving the leading diffusion terms^ ~ b^Y^ and the intercept shift at 0(6); to 
this purpose, expansion (|HT|l is sufficient. We have 



(82) 



'^In principle all diffusion correction terms can be derived using this method. 
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with 



(0)1 



(83a) 

(7) 
(83b) 

(7) 
(83c) 



where, inside the integrals, we have used the same notation for the kernels and their 7-space 
eigenvalues. We restrict our attention to t = to and perform partial integrations to obtain 

G^\toM=j ^ [(Ao + iA;)GW^ + A2Q(GW^)"-(G(°)' 



(i7 
2Tii 



Ao + iAi + iA'2' , 1 A2xi°^" 



+ 



(84) 



where the 7-variable dependence is understood in the A's, G's and x's- Up to this order, the 
maximal energy dependence comes from the cubic pole, which yields a ~ ft^agl"^ dependence. 
The double pole yields instead terms ~ ba^Y which provide the 0{b) correction to ujg- By 
noting that 



G'(°)(y;to,io) 



J 



^^-Xw' (7) 
-1 

(0) ' 



and that a squared Jacobian factor occurs in G^'^\ we obtain the 0(6) correction 

1 



Ao(i) + -A;(i) 



J, 



(85) 
(86) 

(87) 



where actually = 0, because Ai(7) is symmetric for 7 <-> 1 — 7. 

The eigenvalue function Ao(7) is found from the definition in Eq. 1)8 la|) by noting that 
the generalized regularized kernel 



1 



2\ A 



1 



has characteristic function Xl'(7 + f ) + ^r^^ ~ 7 + where 



xf'l(7) 



A 



exp{ Axo(7) + 2^ Xo(7) + 0(A^) 



(88) 



(89) 



and the subscript L [R) refers to the projection with left-hand (right-hand) poles. By proper 
expansion in A we obtain (See App. ^ 

Ao(7) + ^A;(7) = -^{[xg]L(7 + f)-'^X^L(7 + |)-2aoX^/(7) + [7-l-7]} , (90) 
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and finally 



CO 



(1) 



-Or 



_[x^]L(i±^)-a;x^i(i±^)-2aoXrL'(^) 
We note that the expressions of the left projections are (See Add. IX)) 

Xol(7) = V'(l) - iPil) , 



,] + ■■■ 



[xg]L(7)=2[xoL(7)P-V''(7) + ^ 



Xcl(7) 



(91) 

(92) 
(93) 
(94) 



and Xiiil), depending on the resummation scheme, is quoted in App. 1^ 

While the ~ 6^y^ terms exponentiate Aujg and provide a further normalization correc- 
tion j2n], the ~ b'^Y^ terms provide the leading diffusion corrections and occur in 
Considering Eq. H83c|) for t = and performing partial integrations, we obtain at 0{lP') 



AoGW + I'AiGW 



AoGW)'-AiG(o)'(AiGW 



(95) 



This result contains up to a fifth order pole, which can be reduced to a quartic one by partial 
integration, to yield 



^ ' J 27Ti 27Ti 4 _ , 



(0)1 4 



xS'^"(7) Aofffc2xi°^(7,as(fco)) 



27ri 
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[1-5^X^(7)] 



(96) 



The last factor provides the leading diffusion exponent we were looking for. Note that the 
Jacobian factor has been reabsorbed in the fg-derivative of uJs and in the curvature of 
the effective characteristic function: Xeff(^) = x"(l)L-o x s' "^^^^ particular form for the 
generalization of the LO Y^ diffusion term is quite natural when one considers the phys- 
ical mechanism at play: diffusion causes a symmetric spread over a logarithmic range of 
transverse scales of order y/oisX'cs^ ■ The exponent of the evolution at a scale t' is given 

by ujf'\to) + {t' — to)dtoUjf'\tQ). In a first-order expansion of the evolution there is a can- 
cellation between components above and below tQ. But in a second order expansion of the 
evolution, there are corrections from above and below Iq that enter with the same sign, 
~ [±^yasXesYdtoJs'\to)Y]'^. This is precisely the form of (|96l) . 

The analytical treatment given above has its counterpart in the numerical extraction 
of the running-coupling diffusion coefficients presented in |2S|. We illustrate here that the 
method can also be applied to a more general case with an w-dependent resummed NLL 
BFKL kernel. 
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Formally we write the logarithm of the Green's function as a power series in b: 



log G(y; t; to) = Y,V [log G{Y- 1; to)]i , (97) 
i=0 

where the expansion is defined such that as(io) (or optionally some other scale) is kept 
independent of b. We can then write the effective exponent 

ujes (Y ; to) = ^ log G{Y; to- to), (98) 

also as a series in b: 

A 
dY 



^cs{Y; to) = J2 ^'^eff,^ = Yl ^'>^^' ' ^^^^ 



i=0 i=Q 

In practice the power series is determined numerically by carrying out the evolution with a 
generalized 6-dependent coupling af^{k'^), 

4^kfe^) = ,^,,"fl^- (100) 

l + {t- to)basm) 

using several values of 6 = i6b (typically 6b = 0.01 and i ranges from —3 to 3). In the formal 
limit of small 6b, the knowledge of log G(y; toj ^o) ^or n values of b allows one to determine 
the power series up to order b'^~^. 

In Fig. 1^1 we test the analytical prediction for the leading diffusion term ~ Y^ as given 
by Eq. H96() . We show on this plot the term C(;efr,2 from expansion H99|) with the subtracted 
dyY^ term calculated for schemes A and B with scale as{q^) as a function of rapidity Y . 
We clearly see that after the subtraction there is only a linear dependence left, which signals 
presence of the subleading dyY^ terms The numerical value of the diffusion terms is much 
lower in the resummed models than in the LL BFKL equation. For example the coefficient of 
the leading ~ Y^ term, see (|96|) in the LL BFKL case is about 8 times larger than the one in 
the resummed models. As a consequence the regime in which the solution is perturbative is 
much broader in the case of the NLL BFKL. One can see this by studying the contour plots 
in Fig. (|14|) . as will be discussed in more detail in the next section. In particular one finds. 
Fig. (fllk ). that the region where the LL solution is insensitive to non-perturbative results 
is much smaller than in Figs. (|14b.c.d) with the resummed evolution. This result is quite 
encouraging as far as the phenomenological predictions for high energy processes with two 
hard scales are concerned. 

In principle, one could extend our procedure to extract the Y"^ terms too, as has been 
done in ref. j22j for the case of the LL BFKL with running coupling. However, the analytical 
calculation here would be quite involved, since these terms originate from a number of different 
sources, i.e. they come both from (|84j) and (|95|) . and moreover they mix with the terms coming 
from the normalization. In practice these Y"^ terms are expected to be rather small and not 
as relevant for phenomenology as the leading Y^ terms. 

We restrict therefore ourselves to showing only the 0{b) shift to ujs given by the analytical 
expression Eq. (|91|) and compared with the numerical calculation, see Fig. 1101 There is clearly 
a perfect agreement between the two methods, exhibiting the leading Oq behavior of uj'i'\ 

Finally, we show in Fig. ^2 our numerical evaluation of the sum of the first two terms of 
(Jeff, Eq. (|M|) . that is w'-^-' + ^w*^^-', as a function of the coupling constant a^. The correction due 
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Figure 9: The difference between the ujcs,2 coefficient as from \99\) and the leading diffusion 
term calculated from equation 196]) for resummed kernel in schemes A (dashed line) and 
B(solid line). The b-expansion has been performed around the fixed coupling value Qs(/cq) = 
0.1. 

to the running of the couphng reduces somewhat the value of the intercept, as compared with 
the fixed couphng case {b = 0) , which is shown in Fig. ^ The plot in Fig. ^2 summarizes 
our present understanding of ujg, because the higher order terms ~ 6^0,3^... are beyond 
our present level of accuracy, and are perhaps not really meaningful, given the complex 
y-dependence of (|77j) . 

Note that we do not compare directly with our earlier results for ojg Oil, because they 
are based on a different definition (the saddle-point of an effective characteristic function), 
which is less directly related to the Green's function. Nevertheless, the present results are 
consistent with previous ones to within NNLL uncertainties. 

4.5 NP uncertainties on Green's function 

It is well appreciated nowadays that, even with two hard scales, the ultra- high energy behavior 
of the BFKL Green's function is entirely determined by non-perturbative physics. It is 
only in an intermediate high-energy regime that one is able to make reliable perturbative 
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Figure 10: Wefr,! from. Eq. / 1 .9 .9)) as a function of the coupling as calculated in schemes 
A(dashed) and B(solid). Lines represent analytical evaluation based on Eq. iyi\) . the points 
correspond to the numerical extraction. 

predictions HnHmUHlIll] . 

Traditionally one estimates non-perturbative uncertainties on BFKL evolution by ex- 
amining the sensitivity to variations of the infrared regularization of the coupling. More 
recently we showed that a purely perturbative answer can be defined in the context of the 
6-expansion |23,, with the highest perturbatively accessible rapidity being determined by the 
breakdown of convergence of this expansion. In this section we shall examine both approaches. 

Let us consider a variety of infrared (IR) regularizations of the coupling. Mostly we shall 
use cutoff regularizations, 

as{q^)^a^^{q^)e{q-~k), (101) 

with three different values oi k. It will also be instructive to examine a 'freezing' regulariza- 
tion, 

a,(g^) = af'^(max(g^ p)) . (102) 

We believe this freezing regularization to be somewhat less physical, since it allows diffusion 
to arbitrarily low scales in the infrared, in contradiction with confinement. However for the 
purposes of our general discussion it will be helpful to have it too at our disposal. 
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Figure 11: Sum of the first two terms from Eq. I199\) as a function of the coupling Os calculated 
in schemes A(dashed) and B(solid). 

In all cases a^^ is the perturbative one-loop coupling with nj = 4, chosen such that 
0:^(9 GeV^) = 0.244, and no cutoff is placed on exchanged gluon virtualities. The complete 
set of IR regularizations is summarized in table ^ together with the resulting Pomeron 
properties, both for LL and resummation scheme B (NLLb) evolution. 

The two main Pomeron features that one may wish to study are its analytical structure 
and the power, of asymptotic growth, both shown in tabled It is well known that with a 
cutoff one expects the Pomeron to be a pole, while for a frozen coupling one expects a branch 
cut, giving a exp(a;py) growth. Though these properties are most easily derived for 

LL BFKL and a coupling that runs as as(A;^), they apply quite generally. 

As regards the wp values, a first point to note concerns the results for LL evolution, which 
with cutoffs on a^, are much smaller than the naive expectation of Q!s(A;^)xo(1/2) — the 

5 /3 

difference stems from large aj (and higher) contributions, originally noticed by Hancock 
and Ross (discussed also in i48 j ) . 

For NLLb evolution the difference between the cutoff and frozen coupling evolutions is less 
dramatic because of the smaller value of the 'raw' oog value (Figs. IT] and [TT|) . As a result the 
uncertainty on the properties of the 'Pomeron' are somewhat reduced. It is interesting to note 
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k (GeV) 


a,(P) 


asymptotic growth 


ujp (LL) 


cjp (NLLb) 


1.00 (cutoff) 


0.39 


exp{ujpY) 


0.44 


0.32 


0.74 (cutoff) 


0.46 


exp(u;py) 


0.49 


0.35 


0.50 (cutoff) 


0.62 


exp(a;py) 


0.58 


0.41 


0.74 (frozen) 


0.46 


exp(u;py) 


1.28 


0.46 



Table 1: Our set of infrared regularizations of the coupling, together with the resulting asymp- 
totic 'Pomeron' behavior and up values for LL with running coupling as{q) and NLL-q evolu- 
tion. 

that these values for the Pomeron intercept are not too different from those found for the hard 
Pomeron in 'two-Pomeron' fits to data in [IHI- It is not clear however to what extent this can 
be considered significant, since on one hand non-perturbative aspects of small-x evolution are 
likely to be extensively modified by the true non-perturbative physics, including saturation 
effects; and on the other hand because the two-Pomeron fits involve rather strong simplifying 
assumptions. 

Having examined the asymptotic properties of the various infrared regularizations, we can 
now move on to examine the IR sensitivity of 'perturbative' Green's functions. The left hand 
plots of Fig.^l((a) and (c) simply have different rapidity ranges) show G{Y, k—e, k+e) for the 
four infrared coupling regularizations of Tab.^ The transverse momentum k = 4.5 GeV is 
chosen lower than in the plots of section in order enhance the sensitivity to the IR region. 
For reference we also include the uncertainty band due to renormalisation scale uncertainty. 
The discussion that follows will concentrate on the NLLb results, however all the plots of 
Fig. El include also LL results, so as to illustrate the dramatically different IR sensitivity 
between LL and NLLb evolution. 

So let us first consider the three cutoff regularizations (NLLb). One sees that up to 
y ~ 30 they give very similar results. Beyond this point, tunneling occurs (for the lowest 
cutoff), and the three curves start to diverge, indicating that according to this prescription 
the Green's function is no longer under perturbative control. 

When instead one examines the curve with an infrared-frozen coupling, one finds a result 
that at first sight appears paradoxical: the Green's function is somewhat lower than with a 
cutoff regularization, over a wide range of Y in which the cutoff regularization looks relatively 
insensitive to NP effects. Naively one might have expected to see little difference until the 
tunneling point. Our understanding of the observed behavior is that it is connected with 
the use of as{q^) in Eq. (|l8|) . which causes the regularization of the coupling to affect, 
among other things, the virtual corrections of the BFKL equation. Having a larger infrared 
coupling increases the size of the (negative) virtual corrections. In situations where the 
Green's function has a substantially negative second derivative (as it does over a wide range 
of Y) there is an incomplete cancellation with the real contributions (of order 1/Q^), which 
means that a larger infrared coupling leads to smaller preasymptotic growth of the Green's 
function.^ This also explains why the curves with a cutoff IR coupling initially evolve more 
slowly for smaller values of k. 

One could also have imagined more sophisticated IR regularization schemes. For example, 
while maintaining an infrared-frozen coupling, one could have placed an IR cutoff on the 

*One cross-check of this understanding comes from the fact that when evolving with a scale as(fc^) in the 
kernel, differences between cutoff and freezing IR regularizations appear only in the asymptotic Y dependence. 
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Figure 12: (a) Green's function calculated with four different infrared regularizations of the 
coupling, with a renormalisation- scale hand (1/2 < < 2) included for reference for the 
k = 0.74 GeV curve; (b) Green's function calculated in the b-expansion, up to and including 
second, third, fourth and fifth orders in b; (c) and (d) are the same as (a) and (h) respectively 
but on a different scale. In all cases e ~ 0.1/c. 
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Figure 13: Ratios of successive coefficients of IP in the b-expansion of log G(Y, k — €,k + e) 
for y = 5. 



exchanged transverse momentum k. We expect that this would give curves whose initial 
evolution is very similar to that of the IR-frozen coupling case, but whose asymptotic NP 
behavior is a pole, as in the cases with a cutoff on the coupling. 

This confusion arising from this wide range of regularization options was in part the 
motivation for introducing the 6-expansion in [23]. The 6-expansion allows one to define a 
perturbative prediction in close analogy with the prescription that is implicitly contained 
in standard fixed-order perturbative predictions. There, one never has to specify any IR 
regularization. Rather, momentum integrals are implicitly carried out over a perturbative 
fixed-order expansion of the coupling, which is well behaved down to zero momentum. Sensi- 
tivity to non-perturbative effects then manifests itself through the appearance of renormalons 
(see for example the review by Beneke [201)) i-e- factorially divergent coefficients in the series 
expansion for one's observable. 

In a small-x resummation, a pure fixed-order expansion would defeat the purpose of 
the resummation in the first place. However it was shown in [23] that one can expand 
logG in powers of the /^-function coefficient b, and that a truncation of the resulting series 
maintains the advantages of small-x resummation, while providing a prescription for defining 
purely 'perturbative' predictions. This is in addition to its usefulness for studying analytical 
properties of the running-coupling dependence of the Green's function, as has already been 
exploited in Sec. 14.41 

Figs. 112b and 1121 1 show the same Green's function as in Figs. 112b and ll2b . but in trunca- 
tions of the 6-expansion ranging from orders 6^ to b^. One sees how all different truncations 
give fairly similar answers at low Y. But at large Y, the presence of the terms in log G involv- 
ing additional factors of b'^a^Y'^ leads to the splaying out of the different truncations, signaling 
the fundamental limit of the 6-expansion. In certain models (e.g. [^Hl) this is associated with 
the appearance of non-analyticity in b. It is to be noted that this large-y breakdown of the 
6-expansion is not of the renormalon type that is expected in normal perturbative series. 

A detailed study of the figure also reveals that even at low Y the expansion is not entirely 
well-behaved. Indeed successive coefficients of the 6-expansion are all of the same sign and 
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grow quite rapidly, in a way that is suggestive of an infrared renormalon. Infrared renor- 
malons are a factorially divergent behavior of the perturbative series whereby the n^^ ^ 1 
order term is proportional to {asb/p)'""n\ (in simple cases). When interpreted in the language 
of asymptotic series, this translates to an uncertainty on the sum of the perturbative series 
of order {A^/Q^. 

To establish whether it is renormalon behavior that we are seeing, in Fig. we show 
ratios of successive coefficients of 6" in the expansion of log G. The fact that, over a significant 
range of n, one sees a large-n behavior consistent^ with (log G)^") /(log G)^""^) ~ cn, implies 
that it is renormalon behavior. Furthermore by examining a second value of Q one can 
establish that c itself is roughly proportional to a^, c ~ l.Oag. However the constant of 
proportionality, corresponding to a value of p = ^ ~ 0.53, is somewhat surprising, because 
it implies power corrections of order (A/Q)^^, i.e. roughly A/Q. Naively one would have 
expected p = 1 (see also [3^). This difference has yet to be understood, though it should 
be kept mind that significant enhancements of naively expected power-suppressed effects are 
known to be possible due to certain classes of resummation effects [S^. It is interesting 
additionally to note that the formally higher-twist non-perturbative effects that we expect 
for splitting functions in Sec. [5] will also turn out to scale roughly as A/Q rather than A^/Q^. 

Regardless of the precise reason for the unexpected scaling, it can be quite straightfor- 
wardly established that the renormalon behavior is directly connected with the use of as{q'^) 
in the LL part of the kernel; i.e. it has the same origin as the preasymptotic effects that arise 
when modifying the IR regularization of the coupling. Fig. 112b . 

These preasymptotic effects are a feature of BFKL evolution that to the best of our 
knowledge have not been observed before. Given that they are strictly connected to the use 
of as{q'^), they are somewhat model dependent. However the motivations for using as{q'^) 
are quite strong. In particular, as we have mentioned above, this is the scale that is explicitly 
suggested by the form of the NLO corrections; furthermore the appearance of the transverse 
momentum of the emitted gluon as the scale of the coupling is a phenomenon that is well- 
motivated in many other contexts of QCD j53j . 

The appearance of significant preasymptotic NP effects complicates somewhat any at- 
tempt to give a compact summary of NP limits in BFKL evolution. In their absence one 
might have parameterised NP effects at a given transverse scale k, by the rapidity at which 
one loses predictability for the Green's function (e.g. |24l23ll22j l. Instead, we examine contour 
plots. Fig. m of 

Ga{Y;k-e,k + e) 



log- 



(103) 



Gb{Y;k-e,k + e) 

where the subscripts a and b indicate the different non-perturbative treatments in the two 
evaluations of the Green's function. Darker shades indicate good agreement between the 
two evaluations, while lighter shades indicate disagreement. Additionally, to guide the eye, 
we have added explicit contours where the (absolute value of the) log of the ratio is equal 
to 0.1, 0.2 and 0.4, which for brevity we shall refer to as the 10%, 20% and 40% contours 
respectively. 

The first plot. Fig. 114b . given for reference, shows results for LL evolution with two 
different IR cutoffs on the coupling (0.5 GeV and 1 GeV). Preasymptotic effects are fairly 

^Except for the last point — indeed while it is the largest values of n that are the hardest to determine 
accurately with our numerical methods, we have not been able to determine with certainty that the value 
obtained for n = 8 is truly unreliable. Accordingly we have chosen to show the point despite our limited 
confidence in it. 
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Q k = 0.5 GeV / ^k=: GeV I I *3 1^ ^ g.S GeV ^ ^k=: GeV 




k [GeV] k [GeV] 

Figure 14: Contour plots showing the sensitivity of G{Y, k — e,k + e) to the choice of non- 
perturbative regularization, obtained by examining the absolute value of the logarithm of the 
ratio of pairs of regularizations. Darker shades indicate insensitivity to the NP regularization, 
and contours have been drawn where the logarithm of the ratio is equal to 0.1, 0.2 and 0.4. 
Plot (a) shows the result for LL evolution (with as{q)) and two cutoff regularizations (k = 
0.5 GeV and k = 1.0 GeV); (b) shows NLL-q evolution with the same pair of cutoffs; (c) shows 
NLL-Q evolution with truncations of the b-expansion at orders b^ and b^; and (d) shows NLL-q 
evolution, comparing a cutoff regularization (k = 0.74 GeV^ with a b-expansion truncation 
(at order b^). 
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irrelevant here, in part because the asymptotic NP contributions set in quite quickly. The 
contours indicate a linear relation between the maximum perturbatively accessible Y value, 
Ytnaxj and log^, as would be expected if this limit is due to tunneling in the Green's function 
with the lower cutoff. From the simplified version of the tunneling formula |211l22j , 

l^tunncl(fc^)--^^^^^, (104) 

we expect that for asymptotically large k, we should see dY,jii^^/ d log k ~ 2/u>p ~ 3.45. In 
practice, the slope that is measured (for k between 10'^ and 10^ GeV) is about 2.7; given 
that the measurement region is not truly asymptotic, the 20% disagreement between the two 
numbers is not unreasonable. 

Fig. 114b uses the same pair of NP regularizations, but with NLLb evolution. A first 
striking difference is the significant region (lower left-hand quadrant) in which there are 
preasymptotic NP effects at the 20% level. This is connected with the preasymptotic effects 
(due to as{q'^)) mentioned earlier in this section. The second important observation is that the 
rapidity where asymptotic non-perturbative effects become important, l^axj is significantly 
larger than for LL. But as before it is roughly consistent with a manifestation of tunneling in 
the Green's function with the lower cutoff:^'' this time the tunneling formula differs slightly 
from that in |2HI22j. because of the presence of the kinematical constraint in the evolution, 
giving 

W..c.(fe^)-^^-^^^^^^^^^. (105) 

- u;s[k^) 

At very large k one would therefore expect a slope dYraa.x/ ^ogk ~ 2(1 + l/wp) ~ 6.9. The 
measured slope (same k range as above) is roughly 6.1. As for LL evolution, these two results 
are not perfectly compatible, but given that the fe-region is not formally asymptotic, the 
disagreement is not unreasonable. 

As is discussed above, using different infrared regularizations is not the only way of 
gauging non-perturbative effects. Fig. I14h shows what happens if instead we consider two 
truncations of the 6-expansion, at orders and 5^. Once again, for smaller values of k there 
are significant preasymptotic NP effects, though the range of k for which they matter is 
more limited. The upper ('asymptotic') limit on Y due to NP uncertainties also behaves 
differently with the 6-expansion. As was shown in j23], the 6-expansion allows one to reach 
rapidities of the order of the fundamental perturbative limit pite2 5,8,26 . l^ax ~ log^ k"^ / h? . 
This different parametric behavior of l^axj though not directly relevant for phenomenological 
parameter ranges, is evident from the large-fc curvature of the contours, and becomes even 
more so when going to yet larger k. 

The plots so far have shown comparisons of pairs of IR regularizations, or pairs of b- 
expansion truncations. However if we look once again at Fig. 1121 we see that the largest 
preasymptotic 'NP' differences are to be seen when comparing an IR cutoff with the b- 
expansion. Accordingly in Fig. I14l i. we show contours for the ratio of Green's functions 
where one is evolved with a central IR cutoff (k = 0.74 GeV) and the other is determined by 
a b^ truncation of the 6-expansion. This is to be considered as a conservative estimate of the 
impact of non-perturbative effects. 

^°The linear dependence of Fmax on log k only becomes convincingly evident at very large k\ we have limited 
the scale to only moderately large k in order to maintain the visibility of the phenomenologically relevant 
region of k. 
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In this comparison, preasymptotic NP effects are so important at lower k values (below a 
few GeV), that one loses the ability to distinguish them clearly from asymptotic NP effects 
associated with tunneling or diffusion. Only for k > 6 GeV is one able to calculate the 
Green's function over a reasonable range of rapidity (at least up to K = 10) with better than 
20% accuracy. One comes to a similar conclusion if one compares the cutoff and frozen IR 
coupling regularizations, as was illustrated in [H^ . 



5 Resummed anomalous dimension and splitting function 

So far, we have investigated the gluon Green's function in the hard Pomeron regime, in 
which the hard scales fc^,/cQ are of the same order, and — by the 6-expansion method — 
we have isolated diffusion and running coupling effects from the non-perturbative Pomeron 
behavior. In the complementary regime /c^ ^ /cq (or k^ ^> /c^), the collinear properties 
become dominant, and the Green's function is characterized by scaling violations and by the 
corresponding anomalous dimensions. The relation to non-perturbative physics changes also, 
because of the validity of the RG factorization property. By arguments based on the double 
7-representation [TBII^ISI] or on truncated models pni21ll3l] we can state that, for t ^ to, 

Guj{k, ko) = J^uj{k)Pui{kQ) + higher twists , (106) 

where J-^j {^uj) is a solution of the homogeneous equation H23|) which is regular for f — > oo 
(to — oo). While the t-dependence, because of its boundary conditions, is expected to be 
perturbatively calculable, the to-dependence is sensitive to the strong-coupling region and 
to non-perturbative physics, but is factorized so that the standard approach of DGLAP 
evolution pi can apply. We are thus entitled to define 

7res(c^,t) = ^ , 9u{t) = 1 1^ e^'Ul) , (107) 
where fuj{'^) represents J^^^ in 7-space. 



5.1 Resummation by a;-expansion 



The analytical form of the resummed eigenfunction was found in on the basis of the 
(j-expansion — summarized in Sec. 12.41 — which provides the solution 



^(7) = exp 



huj 



X(7) ^ d^X^il) EE x.{l) , 



(108) 



in terms of the eigenvalue function Xi^{l) iii Eqs. (|25|) and H26() . Furthermore, in the "semi- 
classical" regime when 6t > 1/w ^ 1, the behavior of J-^{t) can be found from the saddle 
point estimate 

hojt = Xu.ilu.it)) = Xii^^it)) , (109) 



and the solution is then given by 

J^u^it) = k^Tu^ik) ~ 



^/-27TxUlUt)) 



exp 



dr 7a; (r) 



(110) 
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where the function 7a;(t) satisfies the following identity 

lUt)t - 5^^-(7) = / %{T)dT . (Ill) 
The corresponding gluon anomalous dimension is given by jJU] 

(112) 



xUi) 



1 ix^ii) 

7 2xUl) 



Recall, however, that Eq. H112|) is an acceptable approximation only away from the turning 
point 

X'^^{^{0s,t))=0, (113) 

which is a singularity of H112|) with infinite fiuctuations, and defines the exponent uj = u>s{t) 
at anomalous dimension level. 

Therefore, when u approaches u^sit), one can only rely on the 7-representation (|1()7|1 
in order to define the anomalous dimension past the turning point. This was the method 
followed in |llj (with the choice of scheme in Eq. (|22j) ) in order to provide the resummed 
anomalous dimension and its exponent ujc- In the following, we refer to this this calculation 
as the "w-expansion" result. 

5.2 Practical determination of splitting functions 

Here we are more interested in providing the resummed gluon splitting function directly 
in j;-space, by using the resummation scheme defined by the kernel /Coj and by the corre- 
sponding Green's function. Two methods are available to this purpose. One can exploit 
the 7-representation for the t-dependence on the gluon distribution, and define an anoma- 
lous dimension in w-space as given in Eqs. (|1()7|1()<S|) . To obtain a result in x-space, it is 
then necessary to take the inverse Mellin transform of 7i.cs(w,t). However our formalism for 
calculating the Green's function involves a kernel with higher-order terms in and this 
cannot be straightforwardly represented with a 7-representation, so in order to obtain a 
splitting function within the same 'model' as the Green's function we shall need to resort to 
x-space deconvolution directly from the Green's function, using the numerical method pre- 
sented in This involves calculating the Green's function G{y,t,tQ) and a corresponding 
integrated gluon density 

rQ 

xg{x,Q^)= (fkG^'''-'=^^\\ogl/x,k,ko) , (114) 
and then solving numerically the following equation for the effective splitting function P^q^z, Q^), 



j ^P,^{z,Q^)g{^^^,Q^) . (115) 



(ilogQ2 

In the limit of ^ maxj/^Q, A^}, P(,q{z, Q"^) should be independent of the particular choice 
of fco and of regularization of the coupling, modulo higher-twist corrections. That this is 
true in practice is an important verification of factorization, and provides complementarity 
to analytical 'proofs' based on simplified models. 
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Figure 15: Small-z splitting function determined by two complementary numerical methods 
('y -representation and deconvolution) for the test case of the LL+DGLAP model. For refer- 
ence the pure DGLAP splitting function is also shown. 

As a first step it is interesting to check that the two methods for obtaining sphtting func- 
tions are equivalent. We do this for a 'LL+DGLAP' model (which includes the kinematical 
constraint), namely 

xUl) = 2V'(l) - Hi) -Hl-7 + uj)+ ioAiiuj) (- + , (116) 

\7 1 - 7 + tjy 

where the running coupling is evaluated at scale k. Such a model is of interest because it can 
be fully represented in both the 7-representation, since it has no higher-order terms in as, 
and the Green's function approach, since it is straightforwardly expressed in /c-space. It also 
contains some of the typical sources of potential numerical instability (e.g. the 1/(1 — 2;)+ 
term), making it a powerful 'test-case'. 

Fig. El shows that the effective splitting functions obtained with the two methods are 
nearly identical. The difference between them is of the same order as the higher-twist effects 
that come from varying the regularization of the coupling in the deconvolution method (not 
shown). Also plotted is the 1-loop (LO) pure DGLAP splitting function for comparison. We 
note that at large z one sees the standard 1/(1 — z) behavior in all three curves. 

Having established the validity of the deconvolution approach, one can examine the effec- 
tive splitting functions in the context of the full resummed kernel. We restrict our attention 
to scheme B, given that scheme A is not expected to obey the momentum sum rule. Since 
we are determining a purely gluonic splitting function we take = in the subtracted NLL 
kernel, though we keep = 4 in occurrences of the /? function, so as to maintain a realistic 
running of the coupling. Switching to = 4 in the kernel as well has a relatively small 
effect, cf. Fig.EI 



41 





\ 

\ 


\ 1 1 /fivQ/H rv ^ - 




\ 

\ 


\ LL (ofe(q2)) - - - - 
\ NLLb 

VODGLAP 

\ 

\ 




Q = 4.5Gev'\ 


1 


vOfe(Q) = 0.215 \^ 








ST 






N 




\ 

\ 

^ N ; ■ 
^ N /; 


0.1 


1,1,1,1,' 



10"^° 10"^ 10"^ 10'"^ 10"^ 10° 



z 

Figure 16: Small-z resummed splitting function from resummation scheme B, compared to 
the pure 1-loop DGLAP and BFKL splitting functions (the latter with fixed and running 
couplings). 

Fig. El shows the effective spUtting function for Q = 4.5 GeV. It is compared to the 
1-loop DGLAP sphtting function, and to BFKL splitting functions obtained in the pure LL 
approximation with fixed {ag = as(Q^) ~ 0.215) and running [as{q^)) couplings. 

It is perhaps of interest to discuss first the two LL curves. As can be seen from the figure 
(and as has been discussed extensively elsewhere ^HHHElEniEZI ) j running coupling effects 
alone give strong modifications relative to the fixed-coupling LL splitting function. There 
is a taming of the asymptotic behavior: the cut at ojs = 41og2as — 0.60 is converted to 
a series of poles, the leading one being at u)c — 0.25, with the difference ujs — u)c formally 
of order aj gZlllHlIIIl • The running of the coupling also leads to preasymptotic effects, 
in particular it is associated with a dip at moderately small z. Similar features have been 
discussed by other authors as well, though the details differ: in jJHl the running as as(g^) 
is fully implemented only through to NLL order. In jlH6| lli7j the coupling runs as as{k'^) (a 
NLL difference) and furthermore the use of the Airy approximation in the evaluation of the 
expressions analogous to our Eqs. (|1()7I1()8|) means that their results do not quite correspond 
to an exact solution of Eqs. Q and (|115() . 

From the discussion in section 0] for the Green's function, one expects a further strong 
suppression of the asymptotic growth when going from LL to NLLb — for example (for 6(nj = 
4)) tUg goes from 0.60 to 0.27. However because of non-linearities (and the compensation of 
some double counting), the correction to the splitting function from the combination of 
running coupling and NLLb effects is weaker than would be expected from a simple linear 
combination of the two separate effects. Indeed the final running-coupling, NLLb result for 
(jjc with OsiQ'^) = 0.215 is Wc — 0.18. The preasymptotic dip, to which we return below, is 
also modified in the NLLb resummation, becoming somewhat deeper (about 30% of Os) and 
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z 



Figure 17: The small-z resummed splitting function for k = 0.74 GeV and = 1 together 
with renormalisation scale and IR-regularization uncertainties; the inner band is due to the 
variation of k between 0.5 GeV and 1.0 GeV, while the outer band comes from varying the 
renormalisation scale in the range 1/2 < < 2. Also shown are the splitting function 
obtained with the co-expansion (calculated with a /?o corresponding to Uf = A), the LO 
DGLAP Pgg and the known small-x parts of the NNLO DGLAP Pgg. 

moving to smaller z (~ 10^"^). 

Other important characteristics of the splitting function extracted in scheme B are the 
large-z behavior, which coincides with the expected LO DGLAP result and the value of its 
first moment [lo = 1): the scheme has been constructed such that for fixed coupling, the 
effective characteristic function satisfies asXesil = 0,q;s) = 1. At fixed-coupling, duality 
arguments then automatically lead to the splitting function having a zero first moment 
(to within higher- twist corrections), i.e. validity of the momentum sum rule. More generally, 
for running coupling we expect the momentum sum-rule to hold because at w = 1 the kernel 
is free of leading- twist poles. It is therefore interesting to observe that after full inclusion of 
the IR regularized running coupling, and our rather sophisticated deconvolution approach, 
the numerically derived splitting function of Fig. ^Jdoes indeed have a first moment which is 
zero, to within a few parts in 10^. (We have not so far succeeded in establishing the detailed 
origin of this small departure from zero, though it may well be a higher- twist contribution). 

Given the large difference between the original fixed-coupling LL splitting function and 
the running coupling scheme B result, it is important to establish the order of magnitude 
of potential higher-order and non-perturbative uncertainties. This question is addressed in 
figure 1171 where the scheme B splitting function is shown together with two uncertainty 
bands. The inner band is that associated with the variation of the infrared cutoff k between 
0.5 and 1 GeV, indicating a modest non-perturbative uncertainty.^^ 

^^A more conservative NP uncertainty estimate would consider also an IR frozen coupling. Unfortunately 
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The outer band shows the effect of varying in the range 0.5 < < 2 (a range 
commonly used for fully inclusive quantities). This should give an estimate of the importance 
of potential higher-order corrections. One sees that the main features of the splitting function 
are stable, though at small z the uncertainty grows because different renormalisation scales 
lead to slightly different lOc powers. Another way of investigating higher-order uncertainties 
is to consider the w-expansion of — here recalculated with the same nj convention as 
used for scheme B (n/ = 4 in the /3-function and nj = in the rest of the kernel) and 
transformed to z-space. We recall that the w-expansion is based on the same assumptions as 
scheme B, namely LL-I-NLL BFKL and the requirement of correct LO DGLAP limits. From 
fig- El one sees that down to z ~ 10~^ it agrees with scheme B to within the renormalisation 
scale uncertainties. Below, the NLLb and tj-expansion curves move further apart, essentially 
because their tOc values (0.18 and 0.20 respectively) differ by more than would be expected 
based on the variation. This suggests that for future phenomenological purposes, in the 
very small-z region one might wish to consider the effects of a larger range of x^ variation. 

An aspect of the splitting function that deserves more comment is the dip at moderately 
small z. A priori one may wonder about its origin and indeed whether it might not be some 
form of artefact of our resummation procedure. To help resolve the issue fig. El also shows 
the known small-z part of the NNL DGLAP splitting function (for Uf = 0): 

zPggiz) = as + Bal log ^ , B = + ^ + ^^ -1.549 . (117) 

One sees that the initial decrease of the scheme-B splitting function corresponds closely to 
the decrease of the pure NNL DGLAP splitting function, associated with the B(Xi.\ogl/ z 
term. At a certain point however small-z resummation effects set in and the scheme-B 
structure function starts to rise, giving the characteristic dip structure. The fact that the 
initial decrease of the full Pgg is correlated with an exactly determined (NLLx) piece of the 
NNL DGLAP splitting function suggests that the dip structure is a true feature of the small-z 
splitting function. This belief is reinforced by the observed robustness of the dip structure 
under renormalisation scale and resummation scheme changes (though the depth of the dip 
is subject to some degree of uncertainty).^^ 

An interesting question concerns the impact of the dip on fits to parton distributions. 
Calculations in a (partially) RGI LL model 43 whose effective splitting function also has a 
dip, suggest that it is not incompatible with the available structure function data. Ref. ^Hl 
mentions work in progress on fits involving a resummed splitting function with a dip (actually 
considerably deeper than ours), but detailed results have yet to be presented. It should of 
course be kept in mind that so far we have only presented results for purely gluonic problems 
— phenomenological studies will additionally require a treatment of the quark sector. 

To close off this section, we examine how certain properties of the effective splitting func- 
tion depend on the coupling a^. Fig. 1181 One quantity of interest is the formal small-z 
exponent, Wc, shown in the left-hand plot, together with uncertainty bands from varying the 
IR regularization and the renormalisation scale. One sees that at small a^, regularization 

this leads to numerical instabilities and we are only able to study the case of a coupling frozen down to some 
moderately low scale (below which it is cutoff). From these studies we deduce that including the results from 
a full IR-frozen coupling would roughly double the size of the NP uncertainty band. 

^■^It is worth noting that the dips observed in figs. ll5l and ll6l for the running-couphng LL+DGLAP and LL 
models have different compared to that of scheme B. This is at least in part because the NLLa; terms of their 
low-order expansions differ substantially from the true NLLa; terms contained in scheme B. 
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Figure 18: (a) the small-z exponent, ujc of the effective BFKL splitting function in resumma- 
tion scheme B, compared to the previous oj-expansion result, and to the Green's function oog 
as determined in section \J^ (NLL-q); (h) the position (z^ii^) of the small-z minimum of the 
splitting function and the point (zctoss) below which the resummed splitting function becomes 
larger than the 1-loop DGLAP Pgg. The inner and outer bands have the same meaning as in 
Fig.\T^ 



uncertainties very quickly become negligible, in accord with their expected higher-twist na- 
ture/^ while renormahsation-scale uncertainties decrease much more slowly with ag- Also 
shown, for comparison, are curves for u)c in the w-expansion (quite similar to scheme B) and 
a reproduction of the results of section for Ug to first order in b (here shown with = 4 
in 6, whereas in section n f was uniformly 0). 

Given the late onset of the small-z power growth, other interesting quantities are the 
position of the dip, Zdip, and the point where the effective splitting function becomes larger 
than the plain 1-loop DGLAP splitting function (always defined with = 1), Across- 

Both 

quantities are shown as a function of Og in the right-hand plot of Fig. 1181 As one would per- 
haps expect, as one decreases Ug, one has to go to progressively smaller values of z before the 
BFKL increase of the splitting function becomes visible. In this plot too we note the contrast 
between regularization uncertainties which vanish rapidly with Q and renormalisation scale 
uncertainties which vary much more slowly with Q. 

^■^ Actually the regularization uncertainties seem to decrease roughly as 1 /Q whereas one would have expected 
a behavior — this fact (cf. also the discussion of renormalons for Green's functions in section 0J has yet 

to be understood. 
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6 Inclusion of impact factors 

For a realistic calculation of a physical cross section, high energy factorization requires that 
one specify also the impact factors characterizing the external probes fSec. 12.1'j) . The impact 
factors are known in the LL approximation for a variety of physical processes |27U28| . and 
also in the NLL approximation for virtual photons and for forward jet production j3L)j . 
However, the corresponding expressions are quite involved and still to be implemented in 
numerical algorithms. Furthermore, their accuracy stops at the first non-trivial order in a^. 

The purpose of this section is to show how the resummed scheme for can be extended 
to the corresponding impact factors h's by incorporating subleading corrections due to 1) 
phase-space and threshold effects and 2) leading log collinear singularities. The inclusion of 
the exact NL impact factor expressions |3r)[ l31j is left to a future investigation. 

6.1 Phase space and threshold effects 

Let's first consider deep inelastic scattering 7*(g) + p ^ hadrons in the high energy regime 
u = 2pq = Q'^/xb ^ Q^- According to the analysis presented in Ref. |22j we can factorize 
the LL contribution to the ^*p cross section in the form 



h and / represent the off-shell 7*g* and g*p cross sections in which the virtual gluon has a 
particular polarization. The G function indicates the threshold condition to be satisfied in 
the multi-Regge kinematics (MRK) u ^ ui, U2 ^ , hy the invariants defined in Fig. 119b . 

This threshold condition expresses the fact that the longitudinal part of the momentum 
transfer is small with respect to its longitudinal part (consistency constraint |41ll42lll5] ^. In 
fact, in a frame where the momenta p and q have no transverse component, one has 




(118) 



q = q' — xp (119a) 
k = —zq + zp + k : q ■ k = = p • k . (119b) 



In the last equation, one has to remember the euclidean nature of k: —k^k^ = zzu + k^ . The 
relations among invariants and Sudakov parameters are given by 




(120b) 
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(120c) 



MRK implies v\ ^ hence Eqs. H120all20c|) can be approximated by 
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a) b) 

Figure 19: Kinematic diagrams for: a) deep inelastic scattering; b) 7*7* cross section. The 
variables correspond to (2 times) the scalar product of the corresponding momenta, e.g., 
vi = 2q- k, U2 = 2p- {-k) = 2qo ■ ko, i>2 = 2qo ■ {-k) etc.. 

Therefore, if were replaced by \k^k^\., the function would represent just a phase space 
threshold. As it stands, it yields the consistency condition that the virtuality of the gluon is 
essentially transverse, that is 

^ \k^,k^'\ > ^ , (122) 

i.e., the condition in Eq. (|118|) . The additional thresholds {q + k)"^ ~ i/^ — + fc^ > and 
{p — fe)^ ~ — A;^ > — ensuring the final state particles to be in the forward light cone 
— are implicitly contained in h and / respectively. According to Eq. (|121aj) . we can rewrite 
Eq. (ITTHl) as 

Q) = r ^ , (0!, ^ (!^, a (123) 

Jq2 Ul \Vl J \V J 

= f ^^h(^,Q,k)nz,k) (124) 
in terms of the unintegrated gluon density 

^(.,fe) = |^/(.2,fe)e(.-^) . (125) 

Eq. H124|) is the well known factorization formula for DIS which we present for later conve- 
nience also as a convolution in the invariant "energy variable" vi (Eq. (|12l'{|) ). 
Taking the Mellin transform w.r.t. v/Q"^ = , yields the simpler structure 

'yZ^^iQ) = / ^ hUQ,k)J'Uk) , (126) 

in terms of Mellin transforms^^ of the original factors. 

^''Here we define using as energy scale for vi, at variance with Eq. @, where we used Q\k\ as energy 
scale for i/i. The difference is a multiplicative factor (Q/|fe|)". 
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We remark that the LL behavior of a'^ ^ is determined by the leading (rightmost) singu- 
larity uj = LOp{as) "-^^ of in the cj-plane, while h^} = ho + 0{ui) contributes at LL level 
only through its zero-moment Hq. This amounts to integrating /i(z^i, Q, k) in z^i regardless of 
the I'l -dependence in and identifying in J^: ui = Q^, i.e., z = xb- This shows that the 
details of the phase space effects, in particular those at threshold — evidently ignored in the 
approximation h^j ^ Hq just mentioned — appear only as a NL contribution. On the other 
hand, they are expected to be important when the total energy u has moderately high values. 
Therefore, the w-dependent formulation of impact factors is suitable to describe subleading 
effects coming from the proper treatment of the phase space. 

This applies also to the double fc-factorization formula describing the high energy 7*7* 
cross section: in the MRK u » 1^1,1^2 ^ ^^,^1,^2 ^ Q^iQo'^^'^o threshold condi- 
tion can be applied to any 2^2 sub diagram of Fig. 

UVl VV2 V\V2 V2V\ VV\V2 V\V2 

'^^fcT' ^ ">^>W 

(127) 

The last inequality (obtained by using v > |fc||fco|) shows that the boundary of phase space 
can be very simply described by the combination z/iz/2/(z^|fc| |fco|)- This suggests that we write 
the high energy cross section for photons of polarization A and B (A, B = T,L) as 

a M M = j— — ^^ (128) 



h^{l^l,Q,k)g { "-,k,ko /l^(i^2,Qo,fco) 

where representing the g*g* off shell cross section integrated in the "invariant mass" 
contains the total energy dependence and is constrained by the last of the threshold 
conditions H127|) . 

Eq. (|128|) is just equivalent to the fc-factorization formula in energy space, because 
the convolution in the energy variables can be diagonalized by means of the following Mellin 
transforms: 

ot^QM = r - (^Ya^^MM (129a) 
K{Q,k)= r 'l^(9E.\ h{i^i,Q,k) (129b) 



Q\k\ ^1 V T^l 

GUk, /Co) = - u-g f MM, k, k,) , u= . (129c) 

Jo u \ u J i^\k\\ko\ 



In fact, by using the equality 



QQo Q\k\ i'iU2 Qol^ol 



U V\ Z^|fc||fco| ^2 

and the thresholds (|122() and ()127|) . we obtain 

d^fe Skn 



(130) 
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The choice of a symmetric energy scale uq = QQo leads naturally to symmetric energy scales 
for the individual factors h^j and Q^j, as one can see in Eqs. ()129() . 

The lesson is that even in the double high energy factorization formula, a single Mellin 
variable uj allows one to treat in a proper way the kinematics of the process, in particular the 
threshold effects. This motivates our choice to use w-dependent impact factors and kernel. 



6.2 Collinear resummation of impact factors 

Additional subleading contributions to a'^*^* not taken into account in the Green's function 
are higher order perturbative corrections to the impact factors. Here we want to analyze the 
additional corrections which are important in the collinear limits Q ^ Qq and Q <C Qo- In 
order to keep the discussion as simple as possible, we analyze only the fixed coupling (6 = 0) 
case. 

In Sec. 12. 41 we have shown that we can replace the original RG improved Green's function 
Q of Eq. ^ with Q of Eq. (|3U() — the latter being defined in term of the modified kernel (|32|) 
— up to NNLL differences. Correspondingly one should define impact factors /I's which 
provide the same cross section in the new scheme. 

We begin by considering Eq. ()131() with LO impact factors and with the effective Green's 
function Q. With fixed coupling, both the impact factor and the Green's function are scale 
invariant, and the cross section can be given the integral representation 

c^^^(Q,Qo) = 7^ iV'^.(7) (132) 



QQo J 2TTi VQ^ 
^-(7)^^/^:^(7)^.(7)^5(1-7) (133) 

where we have introduced the Mellin transforms 

/iU7) ^ y (^^J hUQ,k) (134a) 
GM = J d^koi-j^j GUKko) (134b) 

We shall now compare the collinear behavior of (|133|) to that predicted for the total cross 
section in order to find the NLL corrections in h"^ at collinear level. In 7-space the formulation 
of collinear factorization becomes particularly simple in the fixed coupling case, and can be 
stated as follows: the leading logQ^/Qg contribution (Q » Qo) — corresponding to the 
behavior at 7 ~ — ^ for the Mellin transform — to the photon-photon cross section at order 
a^a" is given by 

AB/ \ 'iT^UiV^ s--^ 7a; 7a; 7a; 7- V u 

E^^^^^^-^^ (135) 

in terms of the one-loop anomalous dimensions 7^^ describing the "probability" of the a — > b 
splitting, and of additional "photon-vertex factors" distinguishing the polarization of the 
corresponding photon: 

KJ = 1, K;'' = (7 + |)(l + 0(^))- (136) 
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If we restrict our analysis to gluon emission only, beside the two qq pairs coupling to the two 
photons, in Eq. H135|) all aj but two are just gluons. A further approximation, valid in the 
high energy limit, is to neglect collinear gluon emissions between two quarks belonging to 
the same loop. In fact, this amounts to neglecting q ^ q splittings which are proportional to 
= + 0{lo). Therefore, ai = a„_i = q and 32 • • • a„_2 = g. 

In order to find the collinear behavior of the resummed impact factors, we have to com- 
pare Eq. (|135|) with the fc-factorization formula The collinear behavior of the RGI 
kernel at 6 = is simply 

xUl)^as{x^o+^Xc)-^ (137) 

7+2 

This determines the collinear behavior of the RGI Green's function 



(138) 



n=0 ' ' 2 

The collinear behavior of the LO impact factors with exact kinematics j55| is 

hl'^'Hl) = hl<^'\l - 7) ~ 2a0V|^7S^^:^^ (139a) 

h^^'Hl) = ht^'Hl - 7) ~ 2a./W^l'J^ 2 ■ (139b) 

T + 2 1 + + 4^^^ 

Using the relation 

^q7 = ^2iV,7Sg , (140) 
as 

the Mellin transform of the cross section with LO impact factors assumes the form 
-^^(7)Lo i.p.fac. - ^:^^°^(7)^~U7)/.f (°)(1 - 7) 



4 T/A qg °° / gg \n gq,smg qiyB 



7 + f 7 + f ^oV7 + fy 7 + f 7 + f 

having decomposed the anomalous dimensions relative to the q — > g splitting in a singular 
(xl/uj and a non-singular part: 



Ca 



-+B{uj) 

UJ 



^gq,.ing^^gq,n... (-^43) 



We can see from Eq. ()141() that the structure of Eq. (|135() is reproduced, but in the q — > g 
splitting we are taking into account only the singular part of the anomalous dimension. 

This is not a surprise, because the LO impact factors are by definition coupled to the 
Green's function via a high energy gluon exchange, i.e., a singular splitting. Surprising 
enough is the fact that, using the effective Green's function, it suffices to use upper impact 
factor at LO only, in order to obtain the correct collinear singularities on its side. The reason 
is that the additional factor 

n=0 
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stemming from the resolvent of IC^ (see Eq. pU() ) provides exactly the non-singular splittings 
needed to build up the collinear corrections of the upper impact factor to all orders, in the 
collinear ordering Q ^ k. 

On the other hand, the full expression (|135|) contains the correction factor 1 + ujB{lo) 
w.r.t. Eq. (|141|) . due to the full q — > g anomalous dimension (|142jl . This factor is attributed 
to the lower impact factor in the collinear region ko ^ Qq, so that we can set, at NLL level, 

^5(1 - 7) = h^('Hl - 7) + - 7) 



1 + a< 



7 + 



+ NNLL , (144) 



where the additional term shows right-hand singularities only in the 1 — 7 variable (i.e., 
3f?(l-7) > 1/2). 

Analyzing the opposite ordering Qq ^ Q — thus the leading right-hand singularities in 
the variable 7 (i.e., 3^(7) > 1/2) — yields the modification of the upper impact factor 



B{uj) 



1-7 + 



+ NNLL , (145) 



which differs from Eq. (|144j) by the replacements A<r^ B and 7^1 — 7. 

In conclusion, the high energy cross section can be factorized in the product of the Green's 
function and impact factors h as follows: 



'y^^'iQM = I ^^h^iQ,k)gUk,ko)hEiQo,ko) • (146) 



Such a factorization formula includes the full one-loop anomalous dimensions to all orders, 
the NLL contributions of the Green's function and the NLL phase space effects. Still missing 
are the running coupling effects and the subleading collinear NLL corrections to the impact 
factors. The former could be easily incorporated in the collinear limit on the basis of a 
straightforward generalization of Eq. (|141j) . The latter can be included from the known 
results |31j on the basis of the change of scheme discussed below. 



6.3 Next-to-leading impact factors in the cu-independent formulation 

Here we wish to relate the w-dependent formulation of fc-factorization and of BFKL evolution 
used so far, to the more conventional next-to-leading log expansion of the cross section. We 
shall see that this involves a redefinition of NLL impact factors which is somewhat ambiguous, 
and considerably complicates their collinear structure. For the sake of simplicity, we shall 
provide the relation in the frozen limit. 

We have already encountered the operator relation of the w-dependent Green's function 
to the BFKL one up to NLL order. According to Eq. 1)391) we have 

~ [1 - as{K^ + K^)]-^ [lj - as{Ko + a^K^ + 0{al)] , (147) 

which differs from the pure BFKL-type expansion by the operator factor Ti = [1 — as{KQ + 
K^)]^^. The latter originates from the w-shift (expanded to first order in uj) and from 
the collinear behavior. It was first introduced in 5^ where it was shown to provide energy- 
independent terms, which compensate the symmetrical scale choice sq = kko, so as to provide 
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the effective energy scale s> = max(A;^, k^) for the Green's function. On the other hand, the 
complete cross section includes the impact factors /I's according to Eqs. (|146|1 and H145|l . and 
should be consistent with the NLL parameterisation 



X {k\[i^-as{Ko + asKi)r'\ko){h''^^\Qo,ko) + ash''^'\Qo,ko)) . 



In order to compare Eq. (|146() with (|147|) and H148|) . we see that the factor Ti. has to be 
incorporated in the impact factors. We can do that, given the eigenvalue function ^(7) = 
7^(1 — 7), by defining some function Ti-Lil) such that 7i{'j) = 'Hl{i)'HlO- ~ l)i ^"^^ by 
assigning factor Ti-Lil) — 7)) to impact factor A (B). This decomposition is not unique, 

however, and each solution for TCl corresponds to a choice of fc-factorization scheme |3] in 
the subtraction of the leading term at NLL level in perturbation theory. Furthermore, the 
(j-dependence should be expanded in the /I's also. In 7-space we have 



hf. 



a;=0 



uj=0 



Xo + ri^=o 



(149) 



where in the last line we have exploited that the factor (cj — a^xo) eliminates the high energy 
part of the cross section. 

On the other hand, at NL level the factorization of the 7i factor is achieved by setting 
Wl(7) = 1 + asH{-f) + Oiaj), with 



1 



1 



(1-7)^ 



1 1 



7) 



(150) 



Therefore, the NL contribution to the impact factor in the w-independent expansion becomes 



uj=0 



Xo(7) + /iS(7). 



(151) 



We see from Eq. (jlSlf) that the first two terms both generate higher order singularities 
which — e.g., for A = T — are of type 1/7^ and 1/(1—7)'', ^^'^ come from the w-derivative and 
from multiplication by the singular H term. In order to extract the dynamically interesting 
correction h^^^^ from a perturbative calculation of h^^^^, one has to subtract both terms from 
by a proper choice of H, corresponding to a proper factorization scheme. 



7 Discussion 

In this paper we have presented a formulation of the resummed small-x equation based on 
the renormalisation group constraints. The equation presented here embodies correctly the 
LL and NLL BFKL kernels as well as LL DGLAP evolution. The new equation is very close 
to the formulation proposed previously the main difference being the treatment of the 
collinear terms, which are here treated as (j-dependent terms of the leading kernel. The 
advantage of the small-x equation proposed here is that it shows simple collinear poles only 
and is defined directly in k and rapidity space thus making it easy to study the full gluon 
Green's function rather than just its high-energy exponents. Therefore, after inclusion of the 
impact factors, it can be used in a straightforward way for phenomenological applications to 
processes with two hard scales. 
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In our numerical analysis we have obtained the solutions to this equation in the case of 
fixed and running coupling, we have studied the energy dependence of the Green's function 
both for comparable scales and in the collinear limit and we have extracted the corresponding 
splitting function. 

The analysis of the Green's function has confirmed the fact [21 1126 j that the hard Pomeron 
exponent C0s{as) parametrizes only a transient rapidity dependence of the gluon density, to 
be modified by non-linear diffusion corrections at perturbative level and — beyond some 
critical rapidity — by the non-perturbative Pomeron behavior. Nevertheless, subleading re- 
summation effects not only decrease and stabilize LOg itself (Fig. , but basically weaken the 
non-perturbative Pomeron and considerably increase the range of validity of the perturbative 
behavior (Fig. I14j) by 10-20 units in rapidity compared to leading log expectations. There- 
fore, we are encouraged to trust the resummed perturbative predictions for next generation 
accelerators P^ . 

We have provided resummed results for the gluon splitting function also. This is a purely 
perturbative quantity, as has been verified from its definition as the logarithmic derivative of 
the gluon density, by checking collinear factorization for Q ^ Qq. Here resummation effects 
stabilize the (oscillating) logs hierarchy, and cause a soft departure from the DGLAP result, 
showing a shallow dip in the moderate-x region, followed by the expected power increase in 
the very small-x region, characterized by the splitting function exponent uJc{as)- 

Let us now comment in a little more detail on some interesting features of our results. For 
the high-energy exponents, this work confirms the picture of Ref. [11] . The resummed Green's 
function exponent uJs{cts) turns out to be numerically similar to ag for relevant values of as 
(Fis. [TT|) . This exponent is closely related to the saddle point singularity discussed in [TT] . 
but cannot really be identified with it, due to the presence of diffusion corrections to the 
exponent with the same rapidity dependence. For this reason, the extraction of lOs requires 
the subtraction of the leading diffusion term, which turns out to be small compared 
to leading log expectations. On the other hand, the splitting function exponent uJcictg) is 
substantially below ojs — by about 0.1 for typical — due to well-known running coupling 
effects. 

In addition, we find here interesting preasymptotic effects in the energy dependence of 
the gluon density at comparable scales. In particular, the growth of the NLLb resummed 
density is delayed up to rapidities of order y ~ 4 for ~ 0.2. It is also worth commenting on 
the expectations for the onset of perturbative non-linear (saturation) effects in the evolution. 
These become relevant when the Green's function is of order l/a^ j56j . For our reference 
scales {kt — 5 GeV and ~ 0.2), this translates to Y of order 15, i.e. close to the kinematic 
limit of LHC. 

A special comment is needed for the splitting function's dip in the moderate-x region. It 
is at most a 30% effect with respect to the DGLAP value for ~ 0.2, spread over several 
orders of magnitude in x. It is therefore a shallow dip, associated with several subleading 
effects (notably the small-z terms of the NNLO DGLAP splitting function), and it signals 
a quite moderate departure from pure LO DGLAP evolution. Considering this result and 
the values of ujc{as) just mentioned, the overall picture is that our resummed predictions 
are much closer to low order results than naively expected. In turn, this may provide an 
explanation for the apparent success of low order evolution to fit HERA data, despite the 
size of the effective coupling as(Q^) log(l/x). 

In order to compare the present results to experimental data, we need to include the 
physical impact factors for two-scale processes, and quark evolution for DIS processes. Both 
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issues are well studied. For quarks we can follow |54llllj and for impact factors we have 
shown here how to incorporate the collinear resummation (in the frozen coupling limit). We 
recall that the NL contribution to the impact factors depends on the formulation of evolution 
kernel (e.g. w-dependent or independent). In our w-dependent approach, the impact factors 
have a simple collinear behavior due to the oj shift of leading poles. The relation to the more 
conventional w-independent calculations |SJ is given by the subtraction procedure outlined 
in Sec. 16.81 

On the whole, we have presented here a unified description of small-x deep inelastic 
processes, applicable to both the structure function regime and to the 7*7* kinematics. Re- 
summed results push the validity of perturbative QCD towards higher energies and give 
perhaps a preliminary explanation of the cross-sections' apparent smoothness in the small- 
X regime despite the occurrence, in their description, of large perturbative coefficients and 
various strong-coupling phenomena. 
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A Computation of eigenvalues and projections 

We give here some details of the calculation of the kernel eigenvalues and Mellin transforms 
which are needed for Sees. 13.21 and 14.41 

Let us start from the computation of the eigenvalues of the kernel 



) -\5\q)^H)^{k,k') , (152) 



possibly multiplied by the w-shifting factor > order to obtain the kernel of Eq. (|88|1 

of the text. Such kernels are closely related to the regularized form of the leading BFKL 
kernel (with or without consistency constraint), by possibly including powers of log due to 
the running coupling. 

Note first that, denoting by ^''^^7) the eigenvalue function of Hx, the shifted kernel 
H\ has eigenvalues 

X?'(7+f)+xg'(7-f), (153) 
where the left (right) projections of 7(7) are defined by 

f /" di fir/) ^ /(7n) 

iLfiin) = ± / — 7 = > , (154) 

J»7'^5R7 27rz 7 - 7' ^ 7 - 7n 

7n < 
In > 1 

with the upper (lower) determination of signs and conditions. Note that the last expression 
holds in the particular case of simple pole singularities in the left (right) 7-plane. In fact. 
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such eigenvalues are found by the Fourier transform 



g- 2 (t-t')e(t _ t') + e" 2 (*'-*)e(t' - 1) 



Mt'-t) 



dT 



di 



+ e 



(7-^-7')t 



(155) 



In this equation, the 7' integral can be displaced to the left (right) of 7 + ^ (7 ~ ^) in the 
first (second) term in the r.h.s.. By then performing the r-integral in its convergence region 
we obtain, by the definition (|154|) . the result (|15,'-{j) . 

We then compute x^'^'Ht) itself by applying the kernel H\ to the test function {k''^)'^~'^. 
By using then known integral 

d'^k' 



1 



r(A)r(7)r(i-7-A) 

vr (fc'2)i-7[(fc _ fc')2]i-A r(l-A)r(l-7)r(7 + A) 
and the representation 



r(^ + e) 

Viz) 



exp 



6V(z) + ^eV'(^) + ^e^^iz) + ■ 



we obtain [xo(7) = 2V'(1) - ^"(7) - ^"(1 - 7)] 



x['k7) = T 



r(A)r(7)r(i-7-A) 
A Vr(i-A)r(i-7)r(7 + A) 

Axo(7) + \>?Xo{l) + \>? (2V'"(1) - <(7) - - 7)) + 



~ — ^ exp 



(156) 
(157) 
(158) 



which proves Eq. H89j) of the text. Note the cancellation of the A = singularity between real 
emission and virtual term. 

By expanding H158|) in A we obtain, at order A'^, the BFKL eigenvalue, at order A the 
eigenvalue of the running coupling kernel of Eq. (|38() . and so on. The corresponding kernels 
with consistency constraint are obtained by shifting the left/right projections of Eqs. H153|l 
and (fTM]) . 

The simplest left /right projections are based on the formulas 



Xol(7) = Xoh(1 - 7) = ^(1) - -^(7) = 1^(1 



1-7 



E 



(l + n)(7 + n) 



1 



1 



7 + n 1 + n 



-x'ol(7)=V''(7) = X; 



n=0 



(7 + ny 



, V''(7) + ^''(l - 7) 



TT 



sin (7r7) 



(159) 



(160) 



For the shifted running coupling kernel we need [xq + XqIl and thus [Xo]i(7)- Since 
Xo = XOL + XOR we obtain 



[xI]l{i) = [xol(7)]' + 2[xol(7)xol(1 - 



(161) 



The last term has simple poles in the l.h.s., so that by Eq. ()154|) it is expressed by a sum 
over residues as 



o\ ( \ (-\ \^ 2 tU I Xol(1 + n)(^ -7) 
2[xol(7)xol(1 - 7)]l = xM + 2 g + ' 



(162) 
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where a subtraction at 7 = ^ has been performed to make the series convergent, and Xol(^) 



ixo(^) = 21og2. 



The series in the r.h.s. (with simple poles) is expressed as a combination of [xol]^ (which 
has simple and double poles) and of ^' (with double poles only), as follows: 



vr 



2[xol(7)xol(1 - = [xol(7)] - ^' {l) + ^ 



2[xol(7)xol(1 - i)]r = [xol(1 - 7)]^ - - 7) + ^ 



(163a) 
(163b) 



Eqs. (|163|1 can be proved by checking residues and values at 7 = ^ of l.h.s. and r.h.s. and 
coincides with Eq. (|93jl of the text. 

Finally, the left/right projections of Xiin) as given in Eq. (|^T|) are computed on the basis 
of Fourier transforms of the type (|155|1 by splitting the r integration into the Joo, 0] ([0, — oo[) 
intervals for the L (R) projection. The result is 



Xil{i) = 7:^" (7) +Xql{i) 



+ n(7) -^'^(7) + 



mi'- 



1 \ 67 vr^ 5 n/ 
V7 ^1-77 36 ~ 12 ~ TsiVc 



(164) 



vr 



1 + 7 



+ :i|-3M(7) + fl + 5^^ 



32 



1 / 1 



+ jC(3) 



1 



4 V7^ (1 - if 



1 /I 



2 V7 1-7 



+-(M(7 + 1) + M(7-1) 



11 
16 



M(7) 



where 



n(7) 



dtt^ 



_lLi2(l)-Li2(t) 



00 , / 



l-t 



E 



ra=0 



^l^'{n + l) 
n + 7 



Mil) 



n=0 
1 



7-2 



(n + 7) 
//1+7 



*'(i)-*'(i)-4l) 



(165) 
(166) 
(167) 



The corresponding expressions for xi in scheme A (Eq. (|6U|)) and in scheme B (Eq. (|62j) ) are 
respectively 



xri^(7) = Xil(7) + C^(0)[xol(7 + |) - Xol(7)] 



Xrf(7)=XiL(7) + C(a;)[l+'^vli(u;)] ^ 



-C(0)- 

7 + f 



(168) 
(169) 
(170) 
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